AD-A141  826 


AFWAL- TR-84-4029 


STRENGTH  OF  BOLTED  JOINTS  IN 
LAMINATED  COMPOSITES 


Fu-Kuo  Chang 
Richard  A.  Scott 
George  S.  Springer 


Department  of  Mechanical  Engineering  and  Applied  Mechanics 
The  University  of  Michigan 
Ann  Arbor,  MI  48109 


March  1984 


Final  Report  for  Period  June  1 983-Deceaber  1983 


DTIC 

ELECTE 
JUNS  884 


SELEC 
JUNS 

Y*  B 


Approved  for  Public  Release:  Distribution  Unlimited 


MATERIALS  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
WRI6HT-PATTERS0N  AFB,  OHIO  45433 


84  06  04  044 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used 
for  any  purpose  other  than  in  connection  with  a  definitely  related  Govern¬ 
ment  procurement  operation,  the  United  States  Government  thereby  Incurs  no 
responsibility  nor  any  obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  h,sve  formulated,  furnished,  or  in  any  way  supplied  the  said  draw¬ 
ings,  specifications,  or  other  data,  is  not  to  be  regarded  by  implication 
or  otherwise  as  In  any  manner  licensing  the  holder  or  any  other  person  or 
corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use,  or 
sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ASP/PA) 
and  is  releasable  to  the  National  Technical  Information  Service  (NTIS).  At 
NTIS,  it  will  be  available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


i.WYTsai,  Project 


S .W.  Tsai  ,  ¥rbjec t~f ngTheerl  thief 
Mechanics  and  Surface  Interaction  Branch 
Nonaettalic  Materials  Division 


FOR  THE  COMMANDER 


ffcwaetfcllic  Materials  division 


“If  your  address  has  changed,  if  you  wish  to  be  removed  fro*  our  Mail¬ 
ing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organisation  please 
notify  AfWAl/MlBM,  W-PAfO,  Chip  45433  to  help  us  maintain  a  current  Mailing 
list. 

Copies  of  this  report  should  not  be  returned  unless  return  is  required  by 
security  considerations,  contractual  obligations,  or  notice  on  a  specific 
docuMtnt. 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  Out  Cnfft) 


■« 

,f 


REPORT  DOCUMENTATION  PAGE 

read  instructions 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  1.  GOVT  ACCESSION  NO. 

AFWAL-TR-84  -4029  4b.A/^/f^ 

J  RECIPIENT’*  catalog  number 

4.  TITLE  (vid  Zvtllllt) 

STRENGTH  OF  BOLTED  JOINTS  IN  LAMINATED 

COMPOSITES 

»■  TYPE  OF  REPORT  A  PERIOD  COVCRCO 

Final  Report 

June,  1983-Ueceuiber  1983 

PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHOR^ 

Fu-kuo  Chang 

Richard  A.  Scott 

George  S.  Springer 

CONTRACT  or  GRANT  MUMSERf«J 

F33615-81-C-5050 

».  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Dep  irtmeut  of  Mechanical  Engineering  and  Applied 
Mechanics,  The  University  of  Michigan, 

Ann  Arbor,  Michigan  48109 

10.  PROGRAM  CLEMENT.  PROJECT.  TAM 
AREA  A  WORK  UNIT  NUMBERS 

FY1457-81-02013 

ii.  controlling  office  name  ano  aooreis 

MatLirials  Laboratory  (AFWAL/MLBM) 

Air  Force  Wright  Aeronautical  Laboratories (AF SC) 
Weigh;- Patterson,  AFB,  OH  45433 

1*.  REPORT  DATE 

Mhrch  1984 

11.  NUMBER  OF  PAGES 

175 

1''.  MONITORING  AGENCY  NAME  A  ADORES  u»  Mlt****t  (mm  CmmIIIoI  Oilif) 

_ 

is.  security  class.  r»tiAJ«  *#»« j 

Unclassified 

is«.  utCL  ajsification/oownGRaoTnu. 
schedule 

*4.  WSTMittuTiON  STATEMENT  (allkit 

Approved  for  public  release,  distribution  unlimited 

•  7.  U'JT  RjtnjTiCMl  STATEMENT  (*t  iA«  Mm,ci  mImM Mm  JN»j*  JO.  If  AINwi  In*  JU*. «,t 

•  >.  SUPPLEMENTARY  MQTE4 

IS  «i£  Y  wSB&i  oh  H»PW  *LM  ii  M|»|I«Y  Mi  *y  St-**  MMlItJ 

Composite  Materials 

Failure  hypothesis 

Joints 

halted  Joints 

V 

»«***•  If  fjhHilEfy  ^ 

A  tsethod  is  presented  tor  predicting  the  failure  strength  and  failure  node  of 
pin-loaded  holes  in  t thorite  reinforced  composite  laminates.  Irha  method 
includes  two  steps.  First,  the  stress  distribution  in  the  laminate  is  calcula¬ 
ted  by  the  use  of  a  finite  element  method.  Second,  the  failure  load  and 
failure  mode  are  predicted  by  scans  of  a  proposed  failure  hypothesis  together 
with  the  Yaaada-Suu  failure  criterion.  A  eospueer  code  was  developed,  which 
can  he  used  to  calculate  the  maximum  load  and  the  soda  of  failure  of  joint* 
involving  laminates  with  different  ply  orientations,  different  material  ~  ~ 

•  JAN  *) 


UKClASSIFtEQ 


tCUlRiYV  QF  Vm.'S  i-tst 


I? 


SECURITY  CL ASSI rlCATION  OT  THIS  PAGCf»7»—  Dsim  tmUfJ) 


properties,  and  different  geometries. 

Tests  were  performed,  measuring  the  rail-shear  strength  and  the  characteristic 
lengths  for  Fiberite  T300/l<)34-C  composites.  Tests  were  also  conducted,  measur¬ 
ing  the  failure  strengths  and  failure  modes  of  Fiberite  T300/1034-C  laminates 
containing  a  pin-loaded  hole  or  two  pin-loaded  holes  in  parallel  or  in  series. 

Comparisons  were  made  between  the  data  and  the  results  of  the  model.  Good 
agreement  was  found  between  the  analytical  and  the  experimental  results. 

Using  the  computer  code,  parametric  studies  were  performed,  illustrating  the 
procedures  which  can  be  used  to  si2e  composites  containing  pin-loaded  holes .N 


MASSIFIED 


UOHUTtf  Cl  ASSIGNATION  0*  TtUt  i'-J«  fMKWB 


FOREWORD 


This  report  was  prepared  by  Fu-Kuo  Chang,  Richard  A.  Scott,  and 
George  S.  Springer,  Department  of  Mechanical  Engineering  and  Applied 
Mechanics.  The  University  of  Michigan  for  the  Mechanics  and  Surface 
Interactions  Branch  (AFWAL/MLBH) ,  Nonmetal  lie  Materials  Division, 
Materials  Laboratory,  Air  Force  Wright  Aeronautical  Laboratories, 
Wright- Patterson  AFB,  Ohio.  The  work  was  performed  under  Contract  Num¬ 
ber  F  336lJ^81-C5050,  Project  number  FY1 457-81 -0201 3. 

This  report  covers  work  accomplished  during  the  period  June  1983- 
Oecember  1983. 


II- 


TABLE  OF  CONTENTS 


SECTION  Page- 

I .  INTRODUCTION  1 

II.  PROBLEM  STATEMENT  3 

III.  STRESS  ANALYSIS  7 

3.1  Governing  Equation  6 

3.2  Boundary  Conditions;  Single  Hole  and 

Two  Holes  in  Parallel  13 

3.3  Boundary  Conditions;  Two  Holes  in  Series  17 

3.4  Finite  Element  Analysis  20 

3.4.1.  Method  of  Solution,  Single  Hole  and 

Two  Holes  in  Parallel  25 

3.4.2.  Method  of  Solution;  Two  Holes 

in  Series  26 

IV.  PREDICTION  OF  FAILURE  34 

4.1  Failure  Criterion  34 

4.2  Failure  Hypothesis---Characteristic  Curve  35 

4.3  Solution  Procedure  36 

V.  NUMERICAL  SOLUTION  41 

VI .  EXPERIMENT  47 

6.1  Measurement  Procedure  for  the  Laminate 

Shear  Strength  S  47 

6.2  Measurement  Procedure  for 

The  Characteristic  Length  Rfc  49 

6.3  Measurement  Procedure  sor 

The  Characteristic  Length  R  50 

6.4  Strength  of  Mechanically  Fastened 

Composite  Joints  53 

6.5  Specimen  Preparation  54 


TABLE  OF  CONTENTS  (Concluded) 


VII.  MEASUREMENT  OF  S,  Rfc ,  AND  Rc  58 

7.1  Rail  Shear  Strength  S  58 

7.2  Characteristic  Length  R  59 

7.3  Characteristic  Length  R^  64 

VIII.  EXPERIMENTAL  VALIDATION  OF  THE  MODEL  66 

IX.  DESIGN  CONSIDERATIONS  79 

9.1  Interaction  Coefficients  80 

9.2  Numerical  Values  of  Interaction 

Coefficients  82 

9.3  Laminates  With  One  or  Two  Holes  88 

9.4  Laminates  With  Multiple  Holes  92 

9.5  Failure  Mode  97 

X.  SUMMARY  AND  CONCLUSIONS  101 

REFERENCES  103 

APPENDICES 

u 

A  The  Transformed  Reduced  Stiffness  Matrix  0, .  107 

B  The  Coordinate  Transformation  Matrix  T. .  ^  109 

C  The  Mnite  Element  Mesh  Generator  J  110 

D  Shape  Function  Used  in  the  Finite  Element  Code  113 

E  listing  of  a  Sample  of  Input-Output  of  the  115 

Computer  Code  114 

F  Summary  of  Data  for  Calculating  S,  R.,  129 

G  Summary  of  Data  for  Loaded  Holes  139 


LIST  OP  ILLUSTRATIONS 


r 


i 

i 

i 

I 

i 


i 


/ 


Figure  page 

1.  Descriptions  of  the  Problem.  Top:  Single  Hole 
Model;  Middle:  Two  Holes  in  Series;  Bottom: 

Two  Holes  in  Parallel.  4 

2.  Illustration  of  the  Three  Basic  Failure  Modes  6 

3.  Elastic  Laminates  with  One  Hole  (Left),  Two 

Holes  in  Parallel  (Middle),  and  Two  Holes  in 
Series  (Right).  9 

4.  Configurations  Used  in  the  Finite  Element 

Calculations.  14 

5.  Grid  Used  in  the  Finite  Element  Analysis  for  a 

Single  Hole.  Right  Hand  Figure  is  an  Enlarged 
View  of  the  Grid  Around  the  Hole.  21 

6.  Grid  Used  in  the  Finite  Element  Analysis  for 
Two  Holes  in  Parallel.  Right  Hand  Figure  is  an 
Enlarged  View  of  the  Grid  Around  One  of 

the  Holes.  22 

7.  Grid  Used  in  the  Finite  Element  Analysis  for 

Two  Holes  in  Series.  Right  Hand  Figure  is  an 
Enlarged  View  of  the  Grid  Around  Hole.  23 

Q,  Illustration  of  the  Local  Coordinate  System  x. 

and  x,  along  the  Contact  Surfaces  27 


9.  Illustration  of  the  Reversal  of  the  Normal 
Stresses  When  the  Assumed  Contact  Angles  fr., 
and  0a,  are  Greater  than  the  Actual  Contact^ 
Angles^ Lef t ) .  No  Stress  Reversal  Occurs  for 

the  Actual  Contact  Angles  and  8L  (Right).  32 

10.  Variation  of  the  contact  Angle  With  the  Width 
Ratio  for  a  Laminate  Containing  Two  Holes 

in  Series.  33 

n.  Description  of  the  Characteristic  Curve.  36 

12.  Location  of  Failure  (c-U  Along  the 

Characteristic  Curve.  40 

13.  Stress  o2  *iQn9  a. -axis  in  an  Isotropic 

Infinite  Plate  Containing  a  Circular  Hole. 
Comparison  of  Present  Results  with  Theoretical 
Results  Given  by  Timoshenko  [36].  Parameters 
Used  in  Numerical  Calculations:  5-2. 37  ksi, 
D-2R-0.3  in,  W/D-U,  E/D-8.  L/0-2B.  43 


vii 


R 


LIST  OF  ILLUSTRATIONS  (Cont’d) 


14.  Stress  o2  Along  the  Xj-axis  in  an  Isotropic 
Plate  of ^Finite  Width  Containing  a  Loaded  Hole. 
Comparison  of  the  Present  Results  With  the 
Theoretical  Results  Given  by  De  Jong  [24]. 
Parameters  Used  in  the  Numerical  Calculations: 

D*0 . 3  in,  W/D-5,  E/D=4 ,  L/D-14.  45 

15.  Stress  o2  Along  the  x,-axis  in  an  Orthotropic 
Finite  Plate  C 0/90 3  Containing  a  Circular  Hole. 
Comparison  of  the  Present  Results  With  the 
Theoreticla  Results  Obtained  by  Nuismer  and 
Whitney  [34].  Parameters  Used  in  the  Numerical 
Calculations:  Material:  Graphite/Epoxy 
T300/5208 ,  E.»21.4x10J  ksi,  E2-1.6x10j  ksi, 
G,,*0.77x10J  ksi ,Ut5=0.29,  <3=2. 3  ksi, 

D«1  in,  W/D*3 ,  E/D=4,  L/D® 1 4  46 

16.  Schematic  of  Rail  Shear  Test  Fixture.  48 

17.  Fixture  Used  in  Testing  Loaded  Holes  (Base  Plate 

Geometry  Given  in  Figure  18  and  Table  2).  51 

18.  Base  Plates  Configurations  iSee  Figure  17) 

Plate  Thickness  1/4  in.  The  Dimensions  G,  D, 

and  E  are  Given  in  Table  2.  55 

19.  Variation  in  Rail  Shear  Strength  with  the 
Volume  Fraction  of  0  Degree  Plies  of  Cross  Ply 

Laminates,  o  Data.  -  Fit  to  Data.  SSQ 

•50%  Volume  Fraction  of  0  Degree  Plies  •  30 

19400  psi  61 

20.  Characteristic  Length  in  Tension  Rt  as  Function 

of  Hole  Diameter,  width  Ratio,  <  idcPly 
Orientation.  62 

21.  Variation  of  Characteristic  Length  in  Tension  Rt 

With  Hole  Diameter.  Data  are  for  the  Laminate  c 
Configurations  Given  in  Figure  2C.  63 

22.  Bearing  Strengths  of  Fiberite  T300/1Q54-C 
Laminates  Containing  a  Single  Loaded  Hole. 
Comparisons  Between  the  Data  and  the  Results  of 
the  Model.  The  Failure  Modes  Calculated  by  the 
Model  are  the  Same  as  Those  of  the  Data  Unless 
indicated  by  a  Letter  iu  Parentheses  nest  to 

the  Data  Point.  70 


viii 


t  *  *  * 


LIST  OF  ILLUSTRATIONS  (Con'd) 


23.  Bearing  Strengths  of  Fiberite  T3Q0/1034-C 
Laminates  Containing  a  Single  Loaded  Hole. 
Comparisons  Between  the  Data  and  the  Results  of 
the  Modle.  The  Failure  Modes  Calculated  by  the 
Modle  are  the  Same  as  Those  of  the  Data  Unless 
Indicated  by  a  Letter  in  Parentheses  next  to 

the  Data  Point.  71 

24.  Bearing  Strengths  of  Fiberite  T300/1034-C 
Laminates  Containing  a  Single  Loaded  Hole. 
Comparisons  Between  the  Data  and  the  Results  of 
the  Model.  The  Failure  Modes  Calculated  by  the 
Model  are  the  Same  as  Those  of  the  Data  Unless 
Indicated  by  a  Letter  in  Parentheses  next  to 

the  Data  Point.  72 

25.  Bearing  Strengths  of  Fiberite  T300/1034-C 
Laminates  Containing  Two  Loaded  Holes  in  Parallel. 
Comparisons  Between  the  Data  and  the  Results  of 
the  Model.  The  Failure  Modes  Calculated  by  the 
Model  are  the  Same  as  Those  of  the  Data  Unless 
Indicated  by  a  Letter  in  Parentheses  next  to 

the  Data  Point.  73 

26.  Bearing  Strengths  vf  Fiberite  T300/1D34-C 
Laminates  Containing  Two  Loaded  Holes  in  Parallel. 
Comparisons  Between  the  Data  and  the  Results  of 
the  Model.  The  Failure  Modes  Calculated  by  the 
the  Model  are  the  Same  as  Those  of  the  Data  Unless 
Indicated  by  a  Letter  in  Parentheses  next  to 

the  Data  Point.  74 

27.  Bearing  Strengths  of  Fiberite  T30Q/1034-C 
Laminates  Containing  Two  Loaded  Holes  in  Series. 
Comparisons  Between  the  data  and  the  Results  of 
the  Model.  The  Failure  Modes  Calculated  by  the 
Model  are  the  Same  as  Those  of  the  Data  Unless 
Indicated  by  a  Letter  in  Parentheses  next  eo 

the  Data  Point.  75 

26.  Bearing  Strengths  of  Fiberite  T30Q/IQ34-C 

Laminates  Containing  Two  Loaded  Holes  in  Series. 
Comparisons  Between  the  Data  and  the  Results  of 
the  Model.  The  Failure  Modes  Calculated  by  the 
Model  are  the  Same  as  Those  of  the  Data  Unless 
Indicated  by  a  Letter  in  Parentheses  next  to 
the  Data  Point.  76 


LIST  OF  ILLUSTRATIONS  (Cont'd) 


29.  Bearing  Strengths  of  Fiberite  T300/1034-C 

Laminates  Containing  Two  Loaded  Holes  in  Series. 
Comparisons  Between  the  Data  and  the  Results  of 
the  Model.  The  Failure  Modes  Calculated  by  the 
Model  are  the  Same  as  Those  of  the  Data.  77 

30.  Interaction  Coefficient  for  Two  Holes  in 

Parallel.  Results  of  the  Model.  84 

31.  Interaction  Coefficient  for  Two  Holes  in  Series. 

Results  of  the  Model.  65 

32.  Edge  Interaction  Coefficient.  Results  of 

the  Model.  86 

33.  Side  Interaction  Coefficient.  Results  of  the 

Model.  87 

34.  Description  of  the  Problem  Used  in  Designing 
Laminates  With  a)  Single  Pin-Loaded  Hole,  b) 

Two  Pin-Loaded  Holes  in  Parallel,  c)  Two  Pin- 
Loaded  Holes  in  Series.  89 

35.  Failure  Load  as  a  Function  of  Edo»  Ratio  for 

Laminates  Containing  a  Single  Pi>. -Loaded  Holes. 
Results  of  the  Model.  90 


36.  Failure  Load  (Top)  and  Failure  Load  Per  Unit 
Weight  (Bottom)  of  Laminates  Containing  a  Single 
Pin-Loaded  Hole  ,  Two  Pin-Loaded  Holes  in 
Parallel  and  Two  Pin-Loaded  Holes  in  Series. 

Results  of  the  Model.  93 

37.  Geometry  of  Single  Row  of  Holes  (Top)  and  Two 

Rows  of  Holes  (Bottom).  94 

38.  Failure  Load  (Top)  and  Failure  Load  Per  Unit 

Weight  (Bottom)  of  Laminates  Containing  One  Row 
(Left)  and  Two  Rows  (Right)  of  Pin-Loaded  Holes. 
Results  of  the  Mo^ci.  99 

39.  Failure  Modes  of  Laminates  Containing  a  Single 

Row  (Left)  and  Two  Rows  (Right)  of  Pin-Loaded 
Holes.  Results  of  the  Kodel.  id 

40.  Geometry  of  an  Element  Used  in  the  Finite  Eleme.*>' 
Calculations;  Left:  Element  in  the  x.-x- 
Coordinate  System.  Right:  Element  (Master 
Element)  in  the  Local  ir-s)  Coordinate  System. 

*:  is  the  Coordinate  of  Mode  a  in  the  i 
Direction,  q.  is  the  Displacement  of  Node 
a  in  the  i  Direct  ion  and  (r  ,s  )  are  the 
Coordinates  of  Node  a  in  th£  r°s  Coordinate 
System.  i«*l,2,  a-i.2,3,  or  4. 


114 


LIST  OF  TABLES 


Table  Page 

1.  Input  Parameters  Required  by  Computer  Code  and 

the  Output  Provided  by  Code.  42 

2.  Dimensions  of  the  Base  Plates  in  Figures  17 

and  18  .  All  Units  in  Inches.  56 

3.  Properties  of  Fiberite  T3QQ/1Q34-C  Graphite- 

Epoxy  Composite  57 

4.  The  Characteristic  Length  in  Compression  R  for 

Fiberite  T300/1Q34-C.  Data  Obtained  for  c 
D»0.25  in,  W*2.0  in  ,  L=7.0  in,  £^1,25  in.  65 

5.  Approximate  Differences  Between  the  Experimental 
(P)  and  Calculated  (P  )  Failure  Loads  of 
Laminates  Containing  a  Single  Loaded  Hole. 

The  Numbers  Indicate  Maximum  Differences  (in 
Percent!  for  the  Indicated  Hole  Diameters  and 
Ply  Orientations  78 

6.  Rail  Shear  Strength  S  of  Cross  Ply  [ 0/90 ] 

laminates.  (All  Length  Units  in  inches)  s  130 

7.  Characteristic  Length  intension  R  . 

Ply  Orientation  (( 0/245/90) 3 )s  L  131 

8.  Characteristic  Length  in  Tension  R  . 

Ply  Orientation  (0/(*45)3/90, Js  5  132 

9.  Characteristic  Length  in  Tension  R  . 

Ply  Orientation  (0/(s45)2/905Js  1  i 33 

10.  Characteristic  Length  in  Tension  R  . 

Ply  Orientation  (0/24S/907)g  c  134 

11.  Characteristic  Length  in  Tension  R  . 

Ply  Orientation  ( (902/s60/i30)2)s  1  135 

12.  Characteristic  Length  in  Tension  R  . 

Ply  Orientation  ((0/90)*)  1  136 

v  9> 

13.  Characteristic  Length  in  Tension  R  . 

Ply  Orientation  ((245). )„  1  *37 

©  © 

14.  Characteristic  Length  in  Compression  Rc  i 38 

15.  Data  and  Calculated  Values  for  Joints 

Containing  a  Single  Hole.  ( (0/s45/90)3]s  140 

16.  Data  and  Calculated  Values  for  Joints 


■i 


LIST  OF  TABLES  (Cont'd) 


Containing  a  Single 

Hole.  [0/(±45)3/903]g 

142 

17. 

Data  and  Calculated 
Containing  a  Single 

Values  for  Joints 

Hole.  [0/(±45)2/9053s 

143 

18. 

Data  and  Calculated 
Containing  a  Single 

Values  for  Joints 

Hole.  [0/±45/90?]g 

144 

19. 

Data  and  Calculated 
Containing  a  Single 

Values  for  Joints 

Hole.  [(902/±60/±30)23s 

145 

20. 

Data  and  Calculated 
Containing  a  Single 

Values  for  Joints 

Hole.  C(0/90)63s 

147 

21. 

Data  and  Calculated 
Containing  a  Single 

Values  for  Joints 

Hole.  [ (±45) 63g 

149 

22. 

Data  and  Calculated  Values  for  Joints 
Containing  Two  Holes  in  Parallel. 

[ ( 0/±45/90) 3  3  s 

151 

23. 

Data  and  Calculated  Values  for  Joints 
Containing  Two  Holes  in  Parallel. 

[ (902/±60/±30)2Js 

152 

• 

CN 

Data  and  Calculated  Values  for  Joints 
Containing  Two  Holes  in  Parallel.  C ( 0/90 ) g 3 s 

153 

25. 

Data  and  Calculated 

Values  for  Joints 

Containing  Two  Holes  in  Parallel.  [(±45)g]g  154 

26.  Data  and  Calculated  Values  for  Joints 
Containing  Two  Holes  in  Series.  [ (0/±45/90)33g  155 

27.  Data  and  Calculated  Values  for  Joints 
Containing  Two  Holes  in  Series. 

t (902/±60/±30)2]s  156 

28.  Data  and  Calculated  Values  for  Joints 

Containing  Two  Holes  in  Series.  [(0/90)g]g  157 

29.  Data  and  Calculated  Values  for  Joints 

Containing  Two  Holes  in  Series.  [(±45)g3s  158 


LIST  OF  SYMBOLS 


B 

D 

E 


Eijkl 

Emn 

e 

eo 

£ 


GH 

GV 

9e 

9h 

9S 

*V 

H 

hP 


Total  Surface  Area  of  Laminate 

Stress  Prescribed  Area 

Stress  Free  Area 

Displacement  Prescribed  Area 

Surface  Along  Symmetric  Axis 

Total  Contact  Surfaces  Inside  Upper  and  Lower 

Holes 

Surface  Area  of  an  Element  g  on  Which  Surface 
Tractions  are  Applied 

Bearing  Stress 

Diameter  of  Hole 

Edge  Distance 

.-r Stic  Moduli 

Reduced  Elastic  Moduli 

Failure  Indicator  (e<1  Non-Failure,  e^l  Failure) 

Maximum  Value  of  e  on  Characteristic  Curve 

Fraction  of  By-Pass  Load  over  Total  Load 

Assembled  Load  Vector 

Distance  Between  Two  Holes  in  Parallel 

Distance  Between  Two  Holes  in  Series 

Edge  Distance  Coefficient 

Parallel  Hole  Interaction  Coefficient 

Side  Interaction  Coefficient 

Series  Hole  Interaction  Coefficient 

Thickness  of  Laminate 

Thickness  of  p-th  Ply 

xiii 


LIST  OF  SYMBOLS  (Cont'd) 


K 


iBka 

igkct 


L 

L 

M 

N 

N 

N 


s 


o 

a 


P 

P, 


r  1 


Stiffness  Matrix  of  g-th  Element 
Assembled  Stiffness  Matrix 
Plate  Length 

Total  Length  of  Steel  Pin 
Number  of  Element 
Number  of  Plies  in  Laminate 
Number  of  Holes  in  a  Row 
Shape  Function 

Unit  Vector  Normal  to  Surface 
Applied  Load 

Load  Carried  by  Pin  (Pins) 

By-pass  Load 

Failure  Load  of  Laminate  Containing  One  Row  of 
Holes 


P  ,  Failure  Load  of  Laminate  Containing  Two  Raw  of 

Holes 

P,  Failure  Load  of  Laminate  (width  w)  Containing 

c  Single  Hole  at  the  Center  of  Laminate 

Pr  Failure  Load  of  Laminate  (width  2W)  Containing 

u  Two  Holes  Separated  by  Distance  6H  (6^  a  w) 

Pw  Failure  Load  of  Laminate  (width  2W)  Containing 

w  Two  Holes  Separated  by  Distance  W 

P«.  Failure  Load  of  Laminate  (width  w)  Containing 

b  Single  Hole  at  Distance  f.  from  the  Edge 

P„  Failure  Load  of  Laminate  (Width  W)  Containing 

1  Two  Holes;  One  Located  at  Distance  E  from  the 

Edge,  the  Other  Located  at  the  Center  of  the 
Laminate 

Pv  Failure  Load  of  Laminate  (width  w)  Containing 

Two  Holes  in  Series 

Pu  Maximum  Failure  Load  of  a  Laminate  Containing 

"  Pin-Loaded  Holes 

P*  Failure  Load  Per  Unit  Weight 


xiv 


%  i 


•  a 


LIST  OF  SYMBOLS  (Cont'd) 


P*  Maximum  Failure  Load  Per  Unit  Weight  of  a 

M  Laminate  Containing  Pin-Loaded  ?»oles 

P*  .  Failure  Load  Per  Unit  Weight  of  a  Laminate 

r  Containing  one  Row  of  Holes 

P*r2  Failure  Load  Per  Unit  Weight  of  a  Laminate 

1  Containing  Two  Rows  of  Holes 

P  Failure  Load  Carried  by  Second  Through 

No_/  Next  to  Last  Pins  in  Laminate  Containing  One  Row 

or  Two  Rows  of  Holes 


P 

side 


s 


Failure  Load  Carried  by  the  Two  Pins  Next  to 
the  Sides  in  a  Laminate  Containing  One  Row  or  Two 
Rows  of  Holes 

The  Distance  Between  the  Side  and  the  Adjacent 
Hole 

Transformed  Reduced  Stiffness  Matrix  of  p-th 
Ply 

Nodal  Displacement 
Radiol  Distance 

Radiol  Distance  to  the  Choracteristic  Curve 
Chatocterist ic  Length  for  Tension 
Characteristic  Length  for  Compression 
Laminate  Shear  Strength 

Laminate  Shear  Strength  of  10/90]  Laminate 
with  50  Percent  Volume  Fraction  or  0  Degree 
Fibers, 

Total  Surface  Area  of  Tvo-Dimensional 
Laminate 

Area  of  Element  g 

Surface  Traction  Component 

Surface  Traction  Component  on  Ajj  Surface 

Surface  Traction  Component  on  A^  Surface 

Normal  Stress  on  Hole  Surface  at  0  •  0 

Displacements 


xv 


LIST  OF  SYMBOLS  (Cortt'd) 


V, 


W 

w 


X 


x3 

y 


o 

e 

e 


f 

u 

L 


pc 


P 


s 


V 


o 


Arbitrary  Displacement  Functions 
Total  Volume  of  Laminate 
Volume  of  Element  g 
Width  of  Plate 

The  Combined  Weight  of  Laminate  and  Pins 

The  Weight  of  the  Laminate 

The  Weight  of  the  Pins 

Ply  Tensile  Strength 

Ply  Compressive  Strength 

Coordinate  Along  Fiber  Direction  in  Each  Ply 

Coordinate  Perpendicular  tothe  Loading 
Direction  in  Laminate  Plane 

Coordinate  Opposed  to  the  Loading  Direction  and 
Perpendicular  to  the  x1  Direction 

Coordinate  Perpendicular  to  the  x1  and  x2  Axes 

Coordinate  Perpendicular  to  the  Fiber  Direction 
in  Each  Ply 

Boundary  Curve  of  Hole  on  Which  Surface 
Traction  is  Applied 

Boundary  of  Element  Along  Contact  Regions 

Boundary  Curve  of  Element  g  on  Which 
Surface  Traction  is  Applied 

Strain  Components  in  x^-Xj  Coordinate  System 

Angle  Measured  Counterclockwise  from  x^-axis 

Angie  at  which  Failure  Occurs 

The  Contact  Angle  on  the  Upper  Hole 

The  Contact  Angle  on  the  Lower  Hole 

Density  of  the  Laminate 

Density  of  the  pin 

Volume  Fraction  of  Plies  with  0  Degree  Fibers 

xvi 


LIST  OF  SYMBOLS  (Concluded) 


o^j  Stress  Components  in  the  x^-Xj  Coordinate  System 

ox,oy,oxy  Stress  Components  in  the  x-y  Coordinate  System 


SECTION  I 


INTRODUCTION 


Among  the  major  advantages  of  laminated  composite 
structures  over  conventional  metal  structures  are  their 
comparatively  high  strength  to  weight  and  stiffness  to 
weight  ratios.  As  a  result,  fiber  reinforced  composite 
materials  have  been  gaining  wide  application  in  aircraft  and 
spacecraft  construction.  These  applications  require  joining 
composites  either  to  composites  or  to  metals.  Most 
commonly,  joints  are  formed  by  using  mechanical  fasteners. 
Therefore,  suitable  methods  must  be  found  to  determine  the 
failure  strengths  and  failure  modes  of  mechanically-fastened 
joints.  A  knowledge  of  the  failure  strength  and  failure 
modes  would  help  in  selecting  the  appropriate  size  joint  in 
a  given  application. 

Owing  to  the  significance  of  the  problem,  several 
investigators  have  developed  analytical  procedures  for 
calculating  the  strength  of  bolted  joints  in  composite 

materials.  Among  the  recent  studies  are  those  of  Waszczak 

and  Cruse  (Reference  1),  Qplinger  and  Gandhi  (References  2  &  3),  Agarwal 
(Reference  4),  Soni  (Reference  5),  Garbo  and  Ogonowski  (Reference  6), 

York,  Wilson,  and  Pipes  (References  7  &  8), and  Col  lings  (9).  The 
results  of  these  investigations  apply  only  to  joints  containing  a 
single  hole,  and,  with  the  exception  of  Agarwal's  method,  none  of  the 
previous  methods  can  predict  the  mode  of  failure.  Furthermore,  as  will  be 

discussed  in  Section  VIII,  the  previous  methods  provide  conservative 
results  and  underestimate  the  failure  strength,  often  by  as  much  as  50 
percent. 
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The  first  objective  of  the  investigation  was, 
therefore,  to  develop  a  method  which  a)  can  be  used  to 
estimate  both  the  failure  strength  and  the  failure  mode  of 
pin-loaded  holes  in  composites,  b)  applies  to  laminates 
containing  either  one  pin-loaded  hole  or  two  pin-loaded 
holes  in  parallel,  or  two  pin-loaded  holes  in  series,  c) 
provide  results  with  better  accuracy  than  the  existing 
analytical  methods  and,  d)  can  be  used  in  the  design  of 
mechanically-fastened  composite  joints.  The  second 
objective  was  to  develop  a  "user  friendly"  computer  code 
which  can  be  used  to  predict  the  failure  strength  and 
failure  mode  of  loaded  holes  (joints)  involving  laminates 
with  different  ply  orientations,  different  material 
properties,  and  different  configurations--  including 
different  hole  sizes,  hole  positions,  and  joint  thicknesses. 
The  third  objective  was  to  generate  data  which  can  be  used 
to  assess  the  accuracies  of  analytical  methods. 

The  analytical  model  and  the  corresponding  numerical 
method  of  solution  are  presented  in  Sections  III -VI .  The 
experimental  apparatus  and  procedures  are  given  in  Section 
VII.  The  data,  and  comparisons  between  the  analytical  and 
experimental  results  are  presented  in  Section  VIII.  The  use 
of  the  model  in  the  design  of  joints  is  described  in  Section 
IS. 
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SECTION  II 
PROBLEM  STATEMENT 

Consider  a  plate  (length  L,  width  W,  thickness  H)  made 
of  N  fiber  reinforced  unidirectional  plies.  The  ply 
orientation  is  arbitrary,  but  must  be  symmetric  with  respect 
to  the  x3=»0  plane  (symmetric  laminate).  Perfect  bonding 
between  each  ply  is  assumed. 

Three  types  of  problems  are  being  analyzed  (Figure  1): 
a)  A  single  hole  of  diameter  D  is  located  along  the 
centerline  of  the  plater  b)  Two  holes  of  diameter  D  are 
located  at  equal  distances  from  the  centerline  of  the  plate 
(two  holes  in  parallel);  c)  Two  holes  of  diameter  D  are 
located  along  the  centerline  of  the  plate  (two  holes  in 
series).  A  rigid  pin,  supported  outside  of  the  plate,  is 
inserted  into  each  hole. 

A  uniform  tensile  load  P  is  applied  to  the  lower  edge  of 
the  plate  and  a  uniform  tensile  load  P2  (referred  to  as  the 
"by-pass-  load)  is  applied  to  the  upper  edge.  These  loads 
are  parallel  to  the  plate  (in-plane  loading)  and  are 
symmetric  with  respect  to  the  centerline.  Hence,  the  loads 
cannot  create  bending  moments  about  either  the  x , ,  x-,,  or  Xj 
axes.  Moreover,  for  symmetric  laminates,  in-plane  loading 
and  bending  effects  are  uncoupled.  Transverse  forces, 

(i.e.,  forces  in  the  direction)  are  not  applied,  and 
transverse  displacement  of  the  laminate  is  not  taken  into 
account.  For  example,  a  washer  on  each  side  of  the 


*  r-2 

ttVttt  f tttf «  f 


P  p 


Figur*  1.  Descriptions  of  the  Problem.  Tops  Single  Hole 
Model;  Middle;  Two  Holes  in  Series;  Bottom; 
Two  Holes  in  Parallel. 
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laminate,  supported  by  a  lightly-tightened  ("finger-tight") 
bolt  in  the  hole,  would  ensure  that  there  is  no  transverse 
displacement,  and  that  the  condition  of  two  dimensionality 
is  satisfied  [10]. 

It  is  desired  to  find  : 

1)  the  maximum  (failure)  load  (PM  )  that  can  be  applied 
before  the  joint  fails,  and 

2)  the  mode  of  failure. 

Point  2  refers  to  the  fact  that,  according  to 
experimental  evidence,  mechanically-fastened  joints  under 
tensile  loads  generally  fail  in  three  basic  modes,  referred 
to  as  tension  mode,  shear-out  mode,  and  bearing  mode.  The 
type  of  damage  resulting  from  each  of  these  modes  is 
illustrated  in  Figure  2.  The  objective,  listed  in  point  2 
above,  is  to  determine  which  of  these  modes  will  most  be 
responsible  for  the  failure. 

The  calculation  proceeds  in  three  steps.  For  a  given 
geometry  and  load  : 

1)  the  stress  and  strain  distributions  around  the  hole  are 
calculated, 

2)  the  maximum  (failure)  load  is  predicted, 

3)  the  mode  of  failure  is  determined. 

The  details  of  these  steps  are  presented  in  Sections  HI  and 


IV 
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SHEAROUT 


BEARING 


Figure  3.  Illustration  o£  the  Three  Basic  Failure  Modes 
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SECTION  III 
STRESS  ANALYSIS 

The  calculation  of  stresses  raises  the  issue  of  whether 

a  two  or  three  dimensional  stress  analysis  is  required.  If 

tests  were  to  show  that  the  stacking  sequence  did  not  affect 
* 

the  failure  strength  and  the  failure  mode,  then  a  two 
dimensional  stress  analysis  would  suffice.  Existing 
experimental  evidence  indicates  that  the  stacking  sequence 
is  important  only  when  a)  the  laminate  is  narrow  (and  edge 
effects  are  not  negligible  [11]),  or  b)  the  laminate  is 
unrestrained  laterally  [12].  However,  even  when  the 
stacking  sequence  affects  the  results,  it  seems  to  affect 
the  failure  strength  by  only  10  percent  to  20  percent 
[n-15].  Furthermore,  the  failure  strength  and  the  failure 
mode  seem  unaffected  by  the  stacking  sequence  when  there  is 
a  slight  lateral  constraint  on  the  laminate,  such  as 
provided  by  lightly  tightened  (finger-tight)  bolts 
[10,16,1?]. 

For  these  reasons,  a  two  dimensional  stress  analysis  was 
chosen  for  the  present  work.  As  will  be  demonstrated  in 
Section  VIII,  this  analysis  provides  a  useful  estimate  of 
the  failure  strength  and  the  failure  mode  of  loaded  holes. 

In  addition  to  being  reasonably  accurate,  the  two 
dimensional  analysis  adopted  here  also  provides  a  simple  and 
inexpensive  means  for  calculating  failure  strengths  and 
failure  modes,  making  it  an  attractive  design  aid. 


3.1)  Governing  Equations 


The  stresses  in  the  laminate  are  calculated  on  the  bases 
of  theory  of  anisotropic  elasticity  and  classical-lamination 
plate  theory.  Accordingly,  in  the  analysis,  planes  are 
taken  to  remain  planes,  the  strain  across  the  thickness  is 
taken  to  be  constant  C  e i  -  *  £ (x 1 .Xj) ] ,  and  only  plane  stresses 
are  considered  ^^*023*033*3}  ♦  Under  these  conditions,  in 
the  absence  of  body  forces,  the  condition  of  force 
equilibrium  can  be  expressed  as  [18] 

(1) 

^©2 j / 3* j * 3^2 2^ 0*2“  ® 

In  index  notation  eq.  ID  becomes 

Osj  0  (2) 

*  j  *  J 

Ojj  is  the  stress  component  in  the  plane  normal  to  the 
X|  axis  and  is  in  the  x^  direction.  The  subscripts  i  and  3 
may  have  the  values  of  1  or  2.  Nov  consider  an  elastic 
laminate  of  volume  VQ  containing  a  single  pin-loaded  hole  or 
two  pin-loaded  holes,  as  shown  in  Figure  3.  Loads  are 
applied  over  the  surface  area  A^.  The  displacements  along 
the  surface  area  AR  are  restricted  in  a  manner  described 
subsequently.  The  surface  area  Ag.  is  free  of  applied 
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The  total  surface  area  is 


A»  A^  +  Ar  «-Ap  (3) 

Let  us  denote  by  any  arbitrary  displacement  inside 
the  body.  is  a  test  function.  The  only  requirement  is 
that  uj  be  continuous  and  differentiable.  In  addition, 
along  the  AR  surface,  the  components  of  u*  normal  to  the 
surface  must  be  zero.  By  multiplying  eq.(2)  by  u^  and  by 
taking  the  volume  integral  of  the  resulting  expression,  we 
obtain 


//A,  «ij,j  °i  av  ■  0 


(1) 


By  employing  the  identity 


"  toijui,,j  '  °ijui,j 


(5) 


and  by  utilising  Gauss1  theorem,  eq.  (4)  may  be  written  as 


•'/»  »lj  "j  “i  **  •  Ulvo  °ij  “i.j 


dV  •>  0 


(6) 


where  n  ^  is  the  unit  vector  normal  to  the  surface.  By 
utilising  eq.(3),  eq.(6)  can  be  expressed  as 


Hk,  0i5nj5id*  *  oS3,’i0i  <>*  *  -f-f*p  "Sj  "j  “i 


dA 


Wv„  °»i  'ui.i 


(7) 


On  the  free  surface  Ap  the  stresses  are  zero.  This 
condition  gives 


°ij  nj  5i  <**  * 0  " 

The  forces  per  unit  area  (called  surface  traction)  at  each 
point  of  the  surface  area  A^  are  [18] 

V  °ij  nj  ( 

Equations  (7)-(9)  yield 


;V‘5* 


°ijniui 


dA  -  J/Jv  dv 


no) 


The  stress  is  related  to  the  displacement  through  the  stress 
'Strain  relationship,  which  for  an  elastic  body  is  [18] 


°ij  *  Eijkl  *kl 


<111 


The  subscripts  k  and  1  may  take  on  the  values  of  1  or  2.  In 
order  to  reduce  the  analysis  from  three  dimensions  to  two 
dimensions,  the  reduced  modulus  is  introduced 


ijkl 


an 


*  n 

*  [(h^/W)  0P, 

pat 


an 


(12) 


vhere  h^  is  the  thickness  of  the  p-th  ply,  and  [§]p  is  the 
transformed  reduced  stiffness  matrix  for  the  P-th  ply  [19] 
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(Appendix  A).  The  subscripts  i,j,k,  and  1  are  related  to  m 
and  n  as  follows 


i-j-1  -*• 

m»  1 

k®l»1  - 

n*  1 

H- 

■ 

U. 

■ 

tv> 

4- 

m»2 

k«l*2  * 

n®2 

i  •*  j 

m»3 

k#l  > 

n»3 

(13) 

Note  that  this  reduced  modulus  is  a  constant  and  is 
independent  of  the  thickness  of  the  laminate.  The  strains 
are  related  to  the  displacements  Uj  by  the  expression 

ekl  *  3uk^3xl  +  9uU3xk  )  04) 

By  combining  eqs  ( 10)  —  (14)  we  obtain 


^VEijklui 


jUk,l 


dv  « 


T.u. 

11 


dA  +  //, 


°ijnjui 


dA  (15) 


Since  the  problem  is  treated  as  two  dimensional,  the 
displacements  and,  consequently  the  strains  are  constant 
across  the  laminate.  Hence  the  stresses,  as  defined  by 
eq.(ll),  are  also  constant  across  the  laminate.  However, 
the  on  axis  stresses  in  each  ply  vary  front  ply-to-ply,  and 
are  given  by 


r  °p» ' 

i 

1  c'  1 

L-  (T)[fi]p  < 

1  '2 
Uj 
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where  the  subscripts  x  and  y  represent  the  directions 
parallel  and  normal  to  the  fibers,  respectively.  The  matrix 
[T]  is  the  coordinate  transformation  matrix  given  in 
Appendix  B. 

3.2).  Boundary  Conditions;  Single  Hole  and  Two  Holes  in 
Parallel 


For  problems  involving  a  single  hole  and  two  holes  in 
parallel,  it  is  assumed  that  a  portion  of  the  surface  of 
each  hole  is  subjected  to  a  surface  traction  T^  (Figure  4). 
The  parameter  T^*  is  related  to  the  applied  load.  The 
spatial  distribution  of  T. *  depends  on  the  magnitude  of  the 
applied  load,  on  the  material  properties,  and  on  the 
geometry  in  a  complex  manner.  It  is  extremely  difficult  to 
determine  the  exact  distribution  of  T-*  inside  the  hole 
[20-22],  To  overcome  this  difficulty,  a  cosine  normal  load 
distribution  was  assumed.  With  this  approximation,  a  force 
balance  in  the  x2direction  gives 

Vf  2 

P  ■  P2  +  H  /  (D/2)TX  cos20  de  (17) 

-v/2  2 

where  T„  is  the  normal  stress  at  the  hole  surface  at  0*0. 
x2 

At  any  arbitrary  angle  0  (-ir/2£  e£ir/2),  the  stress  normal  to 
the  surface  is 


Ti 


* 


TX2"i  0056 


(18) 
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Eq(17)  and  (18)  give 

T-*  *  -4  C((P-P2)/irDH)  nt  cose  (19) 

where  P2  is  the  by-pass  load  which  is  a  fraction  f  of  the 
total  load  P 


P2  -fP 


(20) 


The  values  of  either  P  and  P2  or  P  and  f  must  be  specified. 
Thus,  the  surface  traction  on  AL1  can  be  written  as 


Ti*»  “C  (P(1-f )/irDH)ni  cose 


(21) 


The  surface  traction  on  A^2  is 


Ti**«  (P2/HW)  ni  «  (fP/HW)  n. 


(22) 


For  a  single  hole  C  is  equal  to  1  and,  for  two  holes  in 
parallel  it  is  equal  to  1/2.  The  angle  6  varies  from  -w/2 
to  ir/2  in  each  hole.  The  angle  0  is  in  the  x^Xj  Plane« 
and  is  measured  clockwise  from  the  x2  axis  (FigureD.  For 
isotropic  materials,  the  cosine  normal  load  distribution 
(eq.21)  was  found  to  represent  closely  the  actual  load 
distribution  [23].  Calculations  performed  by  previous 
investigators  also  showed  that,  for  composite  materials,  the 
stress  distribution  inside  the  body  is  insensitive  to  the 
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assumed  load  distribution  [1,  6,  24],  Therefore,  eq.  (2T) 
should  suffice  for  the  purpose  of  the  present  analysis, 
which  is  to  determine  the  overall  strength  of  the  joint. 

Equations( 10) , (20) , (21 )  and  (22)  give 

JJJvo  Ei jlcl^i , juk,l  dV  *  //Ali~C(4p(  1-f )/irDH)n.uicose  dA  ♦ 
//A  (fp/HW)  n.  ui  dA  *  // ArO i j  n j  u.  dA  (23) 

We  recall  that  u^  are  functions  that  can  be  selected 
arbitrarily.  The  unknowns  in  eq.(23)  are  the  displacements 
ufe.  Once  u ^  are  known,  the  stresses  at  every  point  can  be 
calculated  from  eqs  (14)  and  (16). 

Solutions  to  eq.(23)  must  be  obtained  subject  to  the 
following  constraints:  a)  Along  the  symmetry  axis  and  along 
the  lower  edge,  displacements  are  allowed  only  in  the 
direction  tangential  to  the  surface.  These  tangential 
displacements  may  occur  freely  without  any  restraints,  b) 
The  intersection  of  the  symmetry  axis  and  the  lower  edge 
must  not  move  (i.e.,  the  intersection  is  rigidly  fixed). 

The  integral  (eq.  23)  over  the  A^  surface  now  applies  to 
the  surfaces  along  the  symmetry  axis  and  along  the  lower 
edge  (Figure  4).  On  these  surfaces,  the  normal  component  of 
the  displacement  and  the  tangential  component  of  the  surface 
traction  are  zero.  Accordingly,  we  have 

//  °ijnj°i  **  "  0 


(24) 
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Equation  (23)  can  now  be  simplified,  and  becomes 


Eijki5i,jVi  dv  ■ 


-  C(4P( 1-f )/trDH)n^u  ^cos0  dA 


(fP/HW)n.ui  dA 


(25) 


The  method  of  solution  of  eq.(25)  is  described  in 
Section  3.4. 


3.3)  Boundary  Condition;  Two  Holes  in  Series 

For  the  problems  of  two  holes  in  series,  the  fractions  of 
the  load  carried  by  each  pin  are  unknown.  To  analyse  the 
the  problem,  it  is  assumed  that  a  uniform  load  distribution 
is  applied  along  the  lower  edge  of  the  plate,  and  it  is 
further  assumed  that  a  rigid  pin  is  inside  each  hole.  The 
assumption  of  the  rigid  pins  implies  that  the  normal 
displacements  are  sero  along  the  contact  surface  (Figure  4). 
The  extent  of  the  contact  surfaces  are  as  yet  unknown  and 
need  to  be  determined. 

The  uniform  load  distribution  on  the  A^,  surface  is 

T^*»  “ ( P/KW )  ni  (26) 

where  H  and  w  are  the  thickness  and  the  width  of  the  place, 
respectively  (Figure  O. 
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Equations  (15),  (22),  (26)  give 


WEi:ki°i,jVidv 


J7x  ~(P/HW)n.iLdA 

AL1  11 


+  //.  ( f P/HW)n- u-dA 

aL2  1  1 


c 


i 


■27) 


As  before,  can  be  selected  arbitrarily,  but  must 
satisfy  the  displacement  boundary  conditions.  Hence,  the 
unknowns  in  eq.(27)  are  the  displacements  u^.  The  solution 
to  eq. ( 27 )  must  be  obtained  with  the  displacement  u^  subject 
to  the  following  constraints: 

a)  Along  the  symmetry  axis,  displacements  are  allowed  only 
in  the  direction  tangential  to  the  surface  (i.e.,  in  the 
*2  direction).  This  tangential  displacement  nay  occur 
freely. 

b)  The  contacts  between  the  rigid  pins  and  the  surfaces  of 
the  holes  are  assumed  to  be  frictionless  and  are  assumed 
to  take  place  through  arcs  bounded  by  the  angles  8y  and 
0^  (Figure  4).  Along  the  arcs  the  surface  displacements 
can  take  place  only  in  the  direction  tangential  to  the 
surface.  Because  o£  the  assumption  of  frictionless 
contact,  this  displacement  may  occur  freely. 

c)  The  radial  displacements  at  the  intersections  of  the 
symmetry  axis  and  the  upper  edge  of  each  hole  are  zero 
(i.e.,  these  intersections  are  rigidly  fixed  ).  This 
corresponds  to  the  rigid  supporting  pins  being  fixed  in 
space. 
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The  integral  over  the  AR  area  now  applies  to  the 
symmetry  axis  and  to  contact  surfaces.  He  express  AR  as 
the  sum  of  two  surfaces 

V  ars  +arc  (28) 

Ars  is  the  surface  area  along  the  symmetry  axis  and  ARC  is 
the  total  contact  surface  inside  the  upper  and  lower  holes. 
Along  the  symmetry  axis,  the  normal  component  of  the 
displacement  and  the  tangential  component  of  surface 
traction  are  zero.  Accordingly,  we  have 


Equation  (27)  gives 


Wvoeljkl“l,juk,ldv  ‘  f\,  * 


(29) 


(30) 


Solution  to  eq.(30)  require  that  the  contact  area  ARC  (i.e., 
the  contact  angles  8y  and  8^  ,  Figure  4)  be  known.  However, 
the  contact  angles  8y  and  8^  are  as  yet  unknown;  therefore, 
these  angles  must  be  determined  before  solutions  for  u^  can 
be  obtained.  Procedures  for  calculating  8y  and  6L  are 
described  in  Section  3.4.  Note  that  the  procedure  was  also 
performed  for  a  single  hole.  Little  difference  was  found 
between  the  predicted  failure  load  and  the  one  predicted 
using  the  stress  boundary  condition. 
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3.4)  Finite  Element  Analysis 


Solutions  to  eq.(25)  and  (30)  were  obtained  by  a  finite 
element  method.  As  a  first  step  in  the  solution  procedure, 
the  volume  V  is  subdivided  into  M  subdomains  of  volume  v 


M 

'  l  Vg 

g^i  g 


(31) 


Eqs.(25)  and  (30)  may  now  be  written  as 


M 

JJ 7/v  EijklQi,-juk,l  dV  *  I//a  Ti*^i  ^  * 

0-1  g  2  1,2  g=i  *Lg.  1  1 

Ni  M 

Ti  ui  **  *  1^*,o0i:nj“i  M 
=1  L9?  g=i  eg  J  J 


(32) 


Tj*  and  T^**  are  the  surface  tractions  given  by  eqs.(2l)  and 

(22)  for  a  single  hole  and  two  holes  in  parallel,  and  by 

eqs.  (22)  and  (26)  for  two  holes  in  series.  A. „  is  the 

ug 

surface  of  an  element  where  the  surface  traction  is  applied. 
At  any  surface  where  load  is  not  applied,  A^gis  xero.  ACg 
is  the  surface  of  an  element  along  the  contact  surfaces. 

For  problems  involving  a  single  hole  and  two  holes  in 
parallel  the  summation  over  A^  is  sero. 

Advantage  is  now  taken  of  the  assumption  that  the 
strains  t*n,c22*  an<3  the  reduced  modulus  E^,  and  the 

stress  (eq.  n)  are  independent  of  the  thickness.  Thus,  the 
three  dimensional  grid,  consisting  of  M  volume  elements,  may 
be 


Finite  Clement  Analysis  for  - 
he  Hand  Figure  it  an  Enlace* 
Ground  the  Hole. 


rigure  6.  (irid  Used  in  the  finite  Eleswsnt  Analysis  Cor 

Two  Holes  in  Parallel,  Right  Hand  Figure  is  an 
Enlarged  View  oC  the  Grid  Around  One  of 
the  Holes. 


•  • 


Grid  Used  in  the  Finite  Element  Analysis  for 
Two  Holes  in  Series.  Right  Hand  Figure  is  an 
Enlarged  View  of  the  Grid  Around  Hole. 
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replaced  by  a  two  dimensional  grid  consisting  of  M  surface 
elements  of  area  s  (Figures  5-7) 


M 

S  ■  [s  (33) 


Equation  (32)  thus  becomes 


M  M 

I//s  Ei-ikiai  iuk  l  ds  *  I/r  dP  * 

j-i  sg  g-i rLgi  1  1 


9 si  Lg2 


9=1  eg 


°i jn jui  dr 


(34) 


where  r^g  and  rLg  are  segments  of  a  line  which  coincide 
with  the  boundary  of  an  element  g  where  the  load  is  applied 
(Figure  5-7).  rcg  denotes  the  boundary  of  an  element  along 
the  contact  regions  bounded  by  6y  and  6^  (Figure  7). 

Isoparametric  4-node  elements  were  used  in  the 
investigation.  The  mesh  was  generated  using  a  mesh 
generator.  This  mesh  generator  was  designed  to 
automatically  generate  grid  sites  around  the  hole  (or  holes) 
in  a  mannet  which  ensures  accurate  resolution  of  the 
stresses  in  the  vicinity  of  the  holes  (Appendix  C).  Smaller 

i 

grids  were  used  around  the  holes  to  obtain  a  better 
resolution  of  the  stresses.  Utilising  the  symmetry  about 
the  *2  axis,  grids  were  placed  on  one  half  of  the  laminate, 
as  illustrated  in  Figures  5-7.  Grids  consisting  of  306,  612 
and  655  elements  were  used  for  problems  involving  a  single 
hole,  two  holes  in  parallel,  and  two  holes  in  series, 
respectively. 
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3.4.1  Method  of  Solution;  Single  Hole  and  Two  Holes  in 
Parallel 


For  problems  involving  a  single  hole  and  two  holes  in 
parallel,  the  displacements  in  each  element  can  be  expressed 
in  terms  of  the  displacements  of  the  four  nodal  points  [25, 
26] 


u-  ■  N  g  • 
i  a  ^ict 

u  •  “N  a  • 
i  a  4ia 


(35) 


The  subscript  a  designates  the  nodal  points  (a  *1,2,3,  or 
4).  is  the  shape  function  described  in  detail  in 
Appendix  D.  q^a  is  the  displacement  at  the  nodal  point  a  in 
the  i  direction. 

Ue  define  a  stiffness  matrix  for  the  q-th  element  as 


R3iBfcot  “  EijklNa,lNa, j  ds 


(36) 


Kgigka  is  an  eight  by  eight  matrix.  The  subscript  a  may 
take  on  the  values  1,  2,  3,  and  4.  The  nodal  displacements 
4*5'e  independent  of  the  surface  and  line 
integrations. 

Accordingly,  eqs(34),  (35),  and  (36)  yield 


8LK<,i6k°  ,k« 

♦  f.  (fP/KW)  fi.  NadF> 

rLg2  1  B 


(37) 


-C(4P/*DH)n4  H-COse  dr 
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The  nodal  displacements  q.D  are  arbitrary  functions  and 

1  P 

hence  eq.  (37)  can  be  written 


^igka  ^i£ 


(38) 


where  the  global  stiffness  matrix  and  the  load  vector 

F  ■ 0  are  given  by 

*  P 

RiBka  *  [  <39) 

9-1 

M 

FiB  =  H  /_  -(4P/i:DH)ni  N.cose  dr 
B  g-1  bgt  * 

♦  f-  (fP/HW)n,N  dr  (40) 

rL92  1  “ 


The  elements  of  and  the  components  of  the  vector  F.  ^ 
are  known;  hence,  q^  can  be  obtained  from  eq.  (38),  using 
the  Gaussian  elimination  method  [27].  Once  are  known, 
the  displacements  are  calculated  from  eq.  (35). 


3.4.2  Method  of  Solution;  Two  Holes  in  Series 


For  problems  involving  two  holes  in  series,  a  local 
coordinate  system  is  employed  along  the  contact  surfaces. 

The  coordinates  of  this  systemix^  and  x'2)  are  everywhere 
normal  and  tangential  to  the  contact  surfaces  as  illustrated 
in  Figure  8. 


I 

z 


I  llustratii 
and  x'  al 


Loi 
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In  this  coordinate  system  the  component  of  ,  u^  and 
o-j  are  denoted  by  the  symbols,  u^  ,  uj ,  and  o'-  j , 
respectively.  These  parameters  in  the  local  coordinate 
system  are  related  to  the  parameters  in  the  fixed  x^  x2 
coordinate  system  by  the  expressions 


u-  =  A-  u 
1  xm  m 

_  f 

V  Aim  um 

°ijnj  *  Aim°mk"k 


(41) 


The  above  transformations  (eq.40)  are  only  used  for  the 
elements  adjacent  to  the  contact  surfaces.  For  all  other 
elements  these  transformations  are  not  employed,  and  we  have 


ui 


-/ 


(42) 


Therefore,  for  elements  adjacent  to  the  contact  surface,  the 
matrix  [A]  is  : 


CA]  - 


cosa  sina 
•sina  cosa 


(43) 


a  is  the  angle  measured  clockwise  from  the  x^  axis  to  the  x1 
axis  [26]  (Figure  8).  For  any  other  elements  which  are  not 


adjacent  to  the  contact  surfaces,  the  matrix  [A]  is 


[A]  » 


1  0 
0  1 


(44) 


Substitution  of  eqs{41)  and  (42)  into  eq.(34)  gives 


M  /  M 

?//_  A,  E*  ji.iA.  u_  -u'  ids-Tf-  (-P/HW)A.  A._ n'ul  dr 

lm  likl  kn  m, i  n.l  lm  in  m  n 

g=i^  g  g=i  L91 

*  ^rt_J£P/Hw)AimAinnm“ndr  *U r  Ai®Aint’mrnr“ndA  U5) 


g=i  L92 


g=i  eg 


On  the  contact  surfaces,  the  normal  component  of  the 
displacements  and  the  tangential  component  of  the  stress  are 
zero  (in  the  new  coordinate  system  x'  and  Xj).  Accordingly, 
the  line  integral  along  the  contact  surfaces  is  zero.  With 
this  simplicif ication,  eq.(45i  becomes 


M 

vs. 


9=1 

!f, 


A.  J Uiui  .jul  ,  ds 

lm  kn  i ]kl  m, 3  n , 1 


9=1  ug  i  g=1  Lg*. 


_/ 


(46) 


The  displacements  at  the  nodal  point  a  are  now  designated  by 
the  symbol  q'  .  wit n  this  notation,  the  displacements  in 

4  Qt 

each  element  become 
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Ni 


/ 


./ 

ui 


No  ^ia 


(47) 


Nq  is  the  shape  function  given  in  Appendix  D.  The 
calculation  now  proceeds  along  the  line  developed  previously 
for  problems  involving  either  a  single  hole  or  two  holes  in 
parallel  (Section  3.4.1).  The  stiffness  matrix  of  the  g-th 
element  is  defined  as 


K^_,_  -  ff„  A.  AWnE  — ,Na  •  dS  (48 

mjjna  JJsg  13AI  a#l  B » D 

As  before,  the  nodal  displacements  are  independent  of 
the  surface  and  line  integrations.  Thus  eqs(48)-U9)  yield 


M  M 

<4,  4a  ■  ^6^rLg1  -tP/HW)*in*innnNe  dr 


9= 

'hLq2(P2'™)Ki*hinWV 


(49) 


The  nodal  displacements  are  arbitrary  functions.  Thus 
eq.(49)  can  be  written  as 


£ 

afina 

where 

a 

afina 


a  m  f 
^na  aft 


*o6na  *nd  ^m£ 


M 

I  k9 

0-1 


aftna 


are  given  by 


(50) 


m 


(51) 
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f»B  ”  t‘/rr„,  '  lP/HM)AimAin"n  Vr 
L9 1 

*  /rLg2  (£P/HW1  Aia  Ain  nn  Vr> 


(52) 


The  elements  of  Kmgna  and  the  components  of  the  vector  Fmg 
are  known,  provided  that  the  components  of  the  matrix  [A]  in 
eq.(50),  are  known.  Hence,  eq(50)  can  be  solved,  once  the 
contact  angles  have  been  determined.  This  can  be 
accomplished  as  follows. 

Values  of  0jj  and  0^,  0°^  and  0aL,  are  assumed  such  that 
0°^  and  0aL  are  greater  than  u/2.  The  displacements  u^  are 
then  calculated  from  eqs  (41),  (42)  and  (47).  Using 
eqs(ii),  (14),  (4t)  and  (42),  the  normal  stresses  along  the 
contact  surfaces  bounded  by  the  arcs  and  0aL  are  then 

calculated.  For  contact  angles  greater  than  the  actual 

contact  angles  compressive  stresses  become  tensile  (stress 
reserval),  as  illustrated  in  Figure  9.  The  angles  9a0  and 
8®^  are  then  decreased  slightly  (by  one  grid  length,  say), 

and  the  stresses  are  calculated  again.  This  procedure  is 

repeated  until  no  reversal  in  sign  of  the  normal  stresses 
occurs  along  the  arcs,  0  to  0y  and  0  to  0L  (i.e.,  both 
contact  surfaces  are  in  compression)  .  These  values,  and 
0^,  are  taken  to  be  the  contact  angles.  As  an  illustration, 
values  of  the  contact  angles  were  calculated  for  Fiberite 
T300/1034-C  composites  with  different  width  ratios.  The 
variation  in  the  contact  angles  with  the  width  ratios  are 
given  in  Figure  10. 


i  *v‘  } 


Figure  9. 


Illustration  of  the  Reversal  of  the  Normal 
Stresses  when  the  Assumed  Contact  Angles  e* 
and  8  ,  are  Greater  than  the  Actual  Contact0 
Angles  (Left).  No  Stress  Reversal  Occurs  for 
the  Actual  Contact  Angles  8y  and  eL  (Right). 


WIOTH  RATIO.  W/D 


SECTION  IV 


PREDICTION  OF  FAILURE 

In  order  to  determine  the  load  at  which  a  joint  fails 
(failure  load)  and  the  mode  of  failure,  the  conditions  for 
failure  must  be  established.  In  this  investigation,  the 
joint  is  taken  to  have  failed  when  certain  combined  stresses 
have  exceeded  a  prescribed  limit  in  any  of  the  plies  along 
a  chosen  curve(denoted  as  the  characteristic  curve).  The 
combined  stress  limit  is  evaluated  using  the  failure 
criterion  proposed  by  Yamada-Sun  (28). 

4.1)  Failure  Criterion 


Numerous  criteria  for  failure  have  been  proposed  in  the 
past  (29,  30-33].  Although  the  concepts  underlying  the 
different  failure  criteria  may  be  different,  the  results  of 
the  various  criteria  are  generally  quite  similar.  In  this 
investigation,  the  Yamada-Sun  failure  criterion  was  adopted 
[28].  This  criterion  is  based  on  the  assumption  that  just 
prior  to  failure  of  the  laminate  .every  ply'  has  failed  due 
to  cracks  along  the  fibers.  This  criterion  states  that 
failure  occurs  when  the  following  condition  is  met  in  any 
one  of  the  plies 


<o„/2t)2  ♦  (oXy/S)2  •  e2 


f  e  «  i 


(53) 
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no  failure 
failure 
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As  indicated  in  eq.  (53)  failure  occurs  when  e  is  equal  to 
or  greater  than  unity.  In  the  above  equation,  o  and  o 

a  xy 

are  the  longitudinal  and  shear  stresses  in  a  ply, 
respectively  (x  and  y  being  the  coordinates  parallel  and 
normal  to  the  fibers  in  the  ply).  S  is  the  rail  shear 
strength  of  a  symmetric,  cross  ply  laminate  [Q/S0]s,  X  is 
either  the  longitudinal  tensile  strength  or  the  longitu/ ! nal 
compressive  strength  of  a  single  piy.  The  tensile  strength 
is  used  when  the  stress  ox  is  in  tension  (o2£>0).  The 
compressive  strength  (X»XC)  is  used  when  ox  is  compressive 
(ox<0). 

4.2)  Failure  Hypothesis-Character ist ic  Curve 

The  hypothesis  is  proposed  here  that  failure  occurs 
when,  in  any  one  of  tne  plies,  the  combined  stresses  satisfy 
on  appropriately-chosen  failu-e  criterion  at  any  point  on  a 
characteristic  curve.  The  the  ncteristie  ^  urve  '.Figure  U) 
is  specified  by  the  expression 

rc(0)  •  0/2  ♦  Rt  ♦  (Bc  -Rt)  cose  t S* > 

The  angle  0.  measured  deck  wise  from  the  x^-Xj  axis,  may 
range  in  value  from  -e/2  to  t/2.  St  and  Pc  are  referred  to 
as  the  characteristic  lengths  for  tension  and  compression. 


Pique©  U.  Description  of  Characteristic  Curve. 
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These  parameters  must  be  determined  experimentally. 

The  concept  of  the  characteristic  length  in  tension  Rt 
was  introduced  by  Whitney  and  Nuismer  [34-37],  In  recent 
years,  several  investigators  utilised  this  concept  in 
analyzing  the  strength  of  loaded  holes.  However,  different 
investigators  used  different  definitions  of  Rfc,  and  employed 
different  procedures  for  determining  the  value  of  Rfc.  As 
will  be  discussed  in  Sections  VI  and  VI I,  the  method 
proposed  here  for  determining  Rfc  differs  from  that  proposed 
by  previous  investigators  [7,  34,  35].  It  is  also  noted 
that  the  characteristic  length  in  compression  Rc  has  not  yet 
been  employed  in  the  strength  analysis  of  loaded  holes. 

In  this  investigation,  the  characteristic  curve  is  used 
together  with  the  Yamada-Sun  failure  criterion.  Accordingly 
(see  eq.  53),  failure  occurs  when  the  parameter  e  is  equal 
to,  or  is  greater  than  unity  at  any  point  on  the 
characteristic  curve 

No  failure  e  <  1  l  at  r  *  rc  (55) 

Failure  e  S  1  i 

It  is  emphasized  chat  the  above  failure  hypothesis  is 
used  here  in  conjunction  with  the  Yamada-Sun  failure 
criterion  (eq.  53).  However,  the  hypothesis  is  general  and 
is  not  restricted  to  the  Yamada-Sun  criterion.  The 
characteristic  curve  proposed  here  may  be  used  with  any 
other  failure  criterion. 


4.3)  Solution  Procedure 


Whether  or  not  a  joint  fails  under  a  given  condition  is 
determined  as  follows.  For  a  given  load 

a)  The  components  of  strains  of  ej2  an<*  ei2  are 

calculated,  using  the  method  of  solution  described  in 
Section  III. 

b)  The  longitudinal  and  shear  stresses  in  each  ply  are 
calculated  using  eq.(l6'. 

c)  The  parameter  e  is  calculated  (eq.53)  along  the 
characteristic  curve. 

d)  If  e  equals  or  exceeds  the  value  of  unity  (e&U  in  any 
ply  along  the  characteristic  curve,  the  joint  is  taken 
to  hove  failed. 

The  procedure  outlined  above  is  used  to  predict  whether 
or  not  failure  occurs  under  a  given  load.  Due  to  the 
assumption  of  a  linear  stress-strain  relationship,  the 
calculated  stresses  are  linearly  proportional  to  the  applied 
load  P.  This  fact,  together  with  Yamada-Sub  failure 
criterion  (eq.53)  gives 


p  e  { 56 ) 

This  relationship  is  utilized  to  determine  the  maximum 
load  which  can  be  imposed  on  the  joint.  For  a  given 
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load  P,  values  of  e  are  calculated  on  the  characteristic 
curve,  as  discussed  above  (points  a-d) .  Note  that  there  are 
two  characteristic  curves  when  there  are  two  holes.  The 
highest  value  of  e  (eQ)  is  then  determined,  and  the  maximum 
load  is  calculated  by  the  expression 

p«ax  *  p/eo  <57> 

The  calculation  procedure  described  in  the  foregoing 
also  provides  the  location  (angle  9j)  at  which  e  first 
reaches  the  value  of  unity  (e»1)  on  the  characteristic  curve 
(Figure12).  A  knowledge  of  9^  provides  an  estimate  of  the 
mode  of  failure.  When  9j  is  small  (efsQ°),  failure  occurs 
by  the  bearing  mode.  When  0{a45°,  failure  is  due  to 
shearout;  when  0£»90°,  failure  is  caused  by  tension. 

In  summary 


-15°  S  S  15° 

bearing  mode 

30°  %  S  60^' 

shearout  mode 

(56) 

75°  S  S  90° 

tension  m^de 

At  intermediate  values  of  9#,  failure  may  be  caused  by  a 
combination  of  these  modes. 


Figure  12. 


Location  of  Failure  (e»l)  Along  the 
Characteristic  Curve. 


SECTION  V 


NUMERICAL  SOLUTION 


A  "user  friendly"  computer  code  (designated  as  BOLT)  was 
developed  which  is  suitable  for  generating  solutions  to  the 
problem  formulated  in  Sections  III-IV.  The  required  input 
parameters  and  the  output  provided  by  the  code  are 
summarized  in  Table  1.  The  input-output  is  illustrated  by 
the  sample  calculations  included  in  Appendix  E. 

In  order  to  assess  the  accuracy  of  the  numerical  method, 
solutions  were  generated  to  problems  for  which  analytical 
solutions  were  available.  Specifically,  stress 
distributions  were  calculated  in  isotropic  plates  containing 
both  unloaded  (open)  and  loaded  holes,  and  in  orthotropic 
plates  containing  unloaded  holes. 

An  analytical  solution  for  the  stress  distribution  in  an 
infinite  (W  .►  «)  isotropic  plate  containing  an  unloaded  hole 
was  given  by  Timoshenko  [38],  The  stress  distribution  in 
such  a  plate  was  also  calculated  by  the  present  method.  The 
parameters  used  in  the  numerical  calculations  are  given  in 
Figure  13.  A  large  width  (W/D*14)  was  uSed  in  the 
calculation  to  approximate  an  infinite  plate.  The  results 
of  the  present  method  and  the  analytical  solution  of 
Timoshenko  are  compared  in  Figure  13.  There  is  excellent 
agreement  between  the  stresses  calculated  by  the  two 
methods. 
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Table  1.  Input  parameters  required  by  the  computer  code  and 
the  output  provided  by  the  code. 

INPUT  PARAMETERS 

1)  Material  Properties 

a)  Longitudinal  and  transverse  Young’s  moduli;  E1  and  Ej 

b)  Shear  modulus,  G12 

c)  Poisson's  ratio, 

d)  Longitudinal  tensile  and  compressive  ply  strength, 

Xfc  and  Xc . 

e)  Rail  shear  strength  of  a  cross  ply  laminate 
I0/90]s,S50 

f)  Characteristic  lengths,  Rfc  and  Rc 

2)  Geometry 

a)  hole  diameter,  D 

b)  thickness,  H 

c)  width,  W 

d)  length,  L 

e)  edge  distance,  E 

f)  distance  between  two  holes,  G  (for  two  holes  only) 

3)  Ply  orientations 
OUTPUT  PARAMETERS 

1)  Failure  load 

2)  failure  mode 


STRESS,  op/cr 


DISTANCE,  *f/R 


Figure  13.  Stress  o.  Along  x.-\*xis  in  an  Isotropic 

In! inite^Plate  Contiining  a  Circular  Mole. 
Comparison  a!  Present  Results  with  Theoretical 
Results  Given  by  Timoshenko  {38)._  Parameters 
Used  in  Numerical  Calculations:  6*2.37  ksi, 
D*2R"0. 3  in,  W/D»14 .  E/D»8,  L/O-28. 
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The  stresses  in  isotropic  plates  containing  loaded  holes 
were  also  calculated.  Plates  of  infinite  and  finite  width 
were  considered.  Calculations  were  performed  for  the 
parameters  given  in  Figure  14. 

As  shown  in  Figure  14,  the  stresses  calculated  by  the 
present  method  are  in  excellent  agreement  with  De  Jong's 
approximate  solution  [24]. 

The  stress  distribution  in  an  orthotropic  plate  of 
finite  width  containing  an  open  (unloaded)  hole  was  also 
calculated.  The  calculations  were  performed  for  a  plate 
with  the  symmetric  laminate  lay  up  of  [0/90]  .  An 
analytical  solution  for  this  problem  was  provided  previously 
by  Nuismer  and  Whitney  [35],  who  modified  Lekhni tski i ' s 
earlier  solution  [39]  for  an  infinite  plate.  The  results 
given  in  Figure  15  show  excellent  agreement  between  the 
stresses  calculated  by  the  present  method  and  by  the 
analytical  solution. 

The  aforementioned  comparisons  indicate  that  the  present 
method  predicts  the  stress  distribution  around  loaded  and 
unloaded  holes  with  high  accuracy. 


STRESS, 


Figure  14.  Stress  o-  Along  the  x.-axis  in  an  Isotropic 

Plate  of  rinite  width  Containing  a  Loaded  Hole. 
Comparison  of  the  Present  Results  with  the 
Theoretical  Results  Given  by  De  Jong  i24j. 
Parameters  Used  in  the  Numerical  Calculations? 
D-0 . 3  in,  W/D»5,  B/O-4,  L/D-14. 
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Stress  o _  Along  the  x,-axis  in  an  Orthotropic 
Finite  Plate  10/90]  Containing  a  Circular  Hole 
Comparison  of  the  Present  Results  with  the 
Theoreticla  Results  Obtained  by  Nuismer  and 
Whitney  [34].  Parameters  Used  in  the  Numerical 
Calculations:  Material:  Graphite/Epoxy 
T300/5208,  E. -21.4x10  ksi,  E  -1.6x10  ksi, 
G,--0.77xl03iksi,u12-0.29,  5-2.3  ksi, 
oil  in.  W/D-3,  E/D=«,  r,/D-14 


SECTION  VI 


EXPERIMENT 

An  experiment  was  performed  to  measure  the  mechanical 
properties  of  composite  laminates  (with  and  without  holes), 
and  the  failure  strengths  and  failure  modes  of  mechanically— 
fastened  composite  joints. 

The  apparatus  and  procedures  used  in  the  tests  are 
described  in  this  section.  A  brief  description  of  the 
procedure  used  to  fabricate  the  test  specimens  is  also 
given. 

6.1)  Measurement  Procedure  for  the  Laminate  Shear  Strength  S 

Rail  shear  tests  were  performed  to  measure  the  laminate 
shear  strength.  Cross  ply  C 0/90] s  laminates  made  of  either 
20  or  24  plies  were  used  in  the  tests.  Laminates  with 
different  volume  fractions  vQ  of  0  plies  were  tested.  vQ  is 
the  number  of  zero  degree  plies  divided  by  the  total  number 
of  plies. 

The  specimens  ranged  from  8  in  to  7.75  in  in  length  and 
2  in  to  1.5  in  in  width.  These  specimen  dimensions  were 
selected  because  it  was  demonstrated  by  previous 
investigators  that  for  such  specimens,  edge  effects  are 
negligible  [40,41].  The  configurations  of  the  rail-shear 
specimens  are  shown  in  Appendix  £. 
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Figure  16.  Schematic  •>£  Rail  Shear  Test  Fixture 
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eight  3/16  in  diameter  holes  were  drilled  along  two 
sides  of  the  specimens,  as  illustrated  in  Figure  16.  The 
specimens  were  placed  between  a  rail-shear  fixture.  The 
geometry  and  dimension  of  this  fixture  are  given  in  Figure 
16.  The  specimen  was  fastened  to  the  rail-shear  fixture  by 
16  bolts.  The  bolts  were  tightened  to  at  least  80  ft-lbf 
torque. 

The  shear  tests  were  performed  by  placing  the  rail-shear 
fixture  into  a  mechanical  testing  machine  and  by  applying  a 
compressive  load.  The  ultimate  failure  load  of  the  laminate 
was  recorded. 
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6.2)  Measurement  Procedure  for  the  Characteristic 

Length  Rfc 

The  characteristic  length  Rt  was  measured  using 
rectangular  specimens  with  an  open  hole  in  the  center  of  the 
specimen.  Tests  were  performed  with  specimens  having 
different  ply  orientations,  different  hole  sires,  and 
different  dimensions  (Appendix  F).  During  each  test,  the 
specimen  was  subjected  to  a  tensile  load  and  the  ultimate 
load  was  recorded,  in  addition,  (after  failure)  the 
specimens  were  inspected  visually  to  establish  the  mode  of 
failure. 

From  the  measured  tensile  strength  the  value  of  »t  was 
determined  as  follows  s 

At  the  failure  load,  the  stresses  in  the  laminate  were 
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calculated,  using  the  model  described  in  Section  III  for  a 
laminate  containing  and  open  hole.  The  stresses  calculated 
in  each  ply  along  the  9*90  line  were  substituted  into  the 
Yamada-Sun  failure  criterion  (eq.  53).  The  point 
along  this  line  was  found  at  which  the  value  e  became  unity. 
The  distance  between  this  point  and  the  edge  of  the  hole  was 
taken  to  be  Rt.  The  values  of  Rt  thus  measured  are 
presented  in  Section  VII. 

6.3 )Measurement  Procedure  for  the  Characteristic 

Length  Rc 

The  characteristic  length  Rc  was  determined  by  the 
following  method.  A  single  hole  was  drilled  into  the 
specimen.  The  position  and  the  diameter  of  the  hole,  the 
specimen  geometry,  and  the  laminate  configurations  used  in 
the  tests  are  given  in  Appendix  F.  The  specimen  was 
inserted  into  a  fixture  shown  in  Figure  i?.  The  top  part  of 
the  fixture  consisted  of  two  3  in  wide  and  5  in  long  steel 
places  (“main  plates*).  A  1.25  in  diameter  and  3.5  in  long 
rod  was  inserted  between  these  plates.  The  rod  was  fastened 
to  the  main  plates  by  bolts.  A  0.5  in  diameter  hole  was 
drilled  along  the  center  line  of  each  main  plate.  1.5  in 
from  the  bottom  edge.  A  0.S  in  dowel  pin  was  inserted  into 
this  hole. 

The  bottom  part  of  the  fixture  consisted  of  two  3  in 
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wide  and  5  in  long  "base  plates."  These  plates  were 
supported  by  the  dowel  pin.  The  material  to  be  tested  was 
placed  between  the  two  base  plates.  A  second  0.5  in 
diameter  dowel  pin  was  passed  through  the  base  plates  and 
the  laminate. 

A  C  clamp  was  placed  around  the  base  plates  near  the 

lower  dowel  pin  and  tightened  by  hand.  The  purpose  of  this 

clamp  was  to  simulate  the  lateral  force  which  would  be 

provided  by  "finger-tight"  bolts  in  the  hole. 

During  the  tests  the  rod  protruding  from  the  main  plates 

was  inserted  into  the  upper  grips,  and  the  laminate  was 

inserted  into  the  lower  grips  of  a  mechanical  testing 

machine.  A  tensile  load  was  applied  by  the  machine  and  the 

ultimate  tensile  strength  was  recorded. 

From  the  measured  tensile  strength,  the  characteristic 

length  R,  was  determined  .  The  stresses  were  calculated 

using  the  model  described  in  Section  III  for  a  loaded  hole. 

A  value  of  R  was  assumed  and  the  characteristic  curve  was 
c 

constructed  in  the  manner  given  in  Section  IV.  The  value  e 
in  the  Yamada-Sun  failure  criterion  (eq.53)  was  determined 
in  each  ply  along  a  segment  of  the  characteristic  curve, 
ranging  from  e- 1 5  to  8*- 15.  The  procedure  was  repeated  for 
different  assumed  values  of  R  until  (in  any  ply)  the  value 
e«1  was  reached  along  the  characteristic  curve  segment 
(-15S0S+15).  This  value  was  then  taken  to  be  R_.  The 

v 

measured  values  of  R  are  presented  subsequently  (Section 

C 

7.3). 


6.4)  Strength  of  Mechanically  Fastened  Composite  Joints 


The  strengths  of  mechanically  fastened  joints  (loaded 
holes)  were  determined  using  rectangular  specimens.  Either 
a  single  hole  or  two  holes  in  parallel  or  in  series  were 
drilled  in  each  specimen.  The  geometries  of  the  specimens 
and  the  laminate  configurations  used  in  the  tescs  are 
described  in  Appendix  G. 

The  test  was  performed  by  placing  the  laminate  into  the 
fixture  described  previously  and  illustrated  in  Figure  17. 

In  each  test  the  same  main  plate  and  the  dowel  pin  were 
uesd.  The  dimensions  of  the  base  plates  were  different, 
depending  upon  the  specimen  configurations.  The  dimension 
of  the  base  plates  are  given  in  Figure  18  and  Table  2. 

During  the  test,  a  lateral  furce  was  applied  with  one  C 
clao.p  to  simulate  the  lateral  force  that  would  be  provided 
by  "finger-tight1*  bolts  placed  in  the  hole.  The  fixture  was 
inserted  into  a  mechanical  testing  machine.  A  tensile  load 
was  applied  and  the  ultimate  tensile  strength  as  recorded. 
After  the  test,  each  specimen  was  inspected  and  the  mode  of 
failure  was  determined. 


The  laminates  were  constructed  from  Fiberite  T30Q/1034-C 
prepreg  tape.  The  panels  were  cured  in  an  autoclave  [43]. 
The  test  specimens  were  cut  by  a  diamond  saw.  The  holes 
were  drilled  with  solid  carbide  drills  for  hole  diameters 
less  than  one  half  inch  and  by  carbide  tip  drills  for  1/2  in 
diameter  holes.  The  nominal  sizes  of  holes  were  0.125  in, 
0.1875  in,  0.25  in,  and  0.5  in.  The  nominal  size  dowel  pins 
were  the  same.  To  provide  a  close  fit,  each  dowel  pin  was 
dressed  down  by  about  O.GQi  in.  The  properties  of  Fiberite 
T3GQ/1Q34-C  are  listed  in  Table  3. 


vmmmmm 


Figure  18.  Base  Plates  Configurations  (See  Figure  17).  Plate  Thickness  1/4  in 
The  Dimensions  G,  D,  and  E  are  Given  in  Table  2. 
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Table  2 

Dimensions 
and  18.  All 

of  the  base  plates 
units  in  inches. 

shown  in 

Figures 

Single  Hole 

D 

G 

E 

plate 

1 

0.5 

1.5 

plate 

2 

0.25 

1.0 

plate 

3 

0.1875 

0.75 

plate 

4 

0.125 

0.5 

Two  Holes 

in  Parallel 

plate 

1 

0.5 

2.5 

1.5 

plate 

2 

0.5 

1.5 

1.5 

plate 

3 

0.25 

1.25 

1.0 

plate 

4 

0.25 

0.75 

1.0 

Two  Holes 

in  Series 

plate 

1 

0.25 

1.25 

0.75 

plate 

2 

0.25 

0.75 

0.75 

plate 

3 

0.1875 

0.625 

0.5 

plate 

4 

0.1875 

0.375 

0.5 

■ 1  w*  . 
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Table  3  Properties  of  Fiberi.te  T300/1034 
composite 


Longitudinal  Young's  modulus,  21 
Transverse  Young's  modulus,  E2 
Shear  Modulus,  G12 
Poisson' s  Ratio  u12 
Longitudinal  tensile  strength,  Xfc 
Longitudinal  compressive  strength,  Xc 
Rail  shear  strength,  S-S5Q 
Characteristic  length  in  tension,  Rt 
Characteristic  length  in  compression,  Rc 


C  graphite/epoxy 


21300000  psi 
1700000  psi 
897000  psi 
0.3 

251000  psi 
200000  psi 
19400  psi 
0.018  i  ■ 
0.07  in 


SECTION  VII 


MEASUREMENTS  OP  S,  Rfc ,  AND  Rc 

Tests  were  performed  to  determine  the  rail-shear 
strength  S  and  the  characteristic  lengths  Rt  and  Rc  of 
Fiberite  T300/1034-C  composites.  These  data  were  generated 
because  they  are  required  in  the  numerical  calculation  of 
the  failure  strength  and  the  failure  mode  of  loaded  holes. 
The  data  obtained  also  indicate  the  sensitivities  of  S,  Rfc, 
and  R  to  such  parameters  as  specimen  geometry  and  laminate 
configuration. 

The  material  properties  used  in  deducing  S,  R,. ,  and  R 
from  the  measured  data  are  listed  in  Table  3.  In  these 
tables  the  values  of  S,  R  ,  and  R  obtained  in  this 
investigation  are  also  included. 

7. 1 )  Rail  Shear  Strength  5 

Rail  shear  tests  were  performed  with  symmetric  cross-ply 
laminates  (0/90)s  having  different  volume  fractions  of  sero 
degree  plies  and  different  geometries.  The  test  conditions 
and  the  test  results  are  summarised  in  Appendix  F.  During 
some  of  the  tests,  cracks  were  observed  near  the  top  and 
bottom  holes  of  the  rail-shear  fixture.  These  cracks 
resulted  in  a  reduction  of  shear  strength.  Specimens  with 
such  cracks  were  not  used  in  calculating  the  rail-shear 
strengths.  Accordingly,  the  rail-shear  strength  of  cross¬ 
ed 


ply  [0/90]s  specimens  having  50  percent  zero-degree  plies  by 
volume  was  found  to  be  S^q  =19.4  ksi. 

The  rail-shear  strength  depends  on  the  volume  fraction 
of  the  zero  degree  plies  in  the  laminates.  At  volume 
fractions  above  50  percent  the  rail  shear  strength  decreases 
(Figure  19).  At  volume  fractions  above  60  percent  the  rsii 
shear  strength  remains  nearly  constant.  Therefore,  when  the 
volume  fraction  of  zero-degree  plies  is  higher  than  50 
percent,  the  rail-shear  strength  corresponding  to  the 
appropriate  volume  fraction  should  be  used  in  calculating 
the  failure  strength  and  the  failure  mode. 

It  was  observed  that  the  shear  stress  to  shear  strain 
relationship  was  nonlinear.  However,  in  the  present  model, 
this  nonlinearity  was  not  taken  into  account.  The 
assumption  of  linear  stress-strain  relation-  may  result  in 
some  error  in  the  calculated  values  of  the  failure  strength 
(Section  VIII),  especially  for  joints  consisting 
predominantly  of  10/90]  and  t i4 5 3 s  laminates. 

7.2)  Characteristic  length  Rt 

The  characteristic  length  St  was  determined  using 
laminates  with  different  geometries  and  different  ply 
orientations.  The  detailed  results  of  the  measurements  are 
given  in  Appendix  F.  The  data  are  summarised  in  Figures  20 
and  21.  Each  data  point  in  these  figures  is  the  average  of 
four  measurements.  Figure  20  shows  the  variations  in 
with  laminate  lay  up.  In  Figure  21,  all  but  one  set  of  data 


are  presented  in  a  single  plot.  The  data  for  [(±45)] 
laminates  were  excluded  from  Figure  21,  because  with  these 
laminates  failure  occurred  not  by  tension,  but  by  tear-out 
along  the  45  dwgree  fibers. 

The  results  in  Figure  20  show  that  the  value  of  Rt 
depends  on  the  hole  diameter,  the  width  ratio  (W/D),  and  the 
ply  orientation.  The  value  of  Rt  increases  with  increasing 
hole  diameter  (Figures  20  and  21).  As  was  discussed 
previously  (Section  4.2),  different  investigators  use 
different  definitions  of  the  characteristic  length.  It  is 
still  noteworthy  that  an  increase  in  characteristic  length 
with  hole  diameter  was  also  observed  by  Whitney  and  Nuismer 
[34,  35]  and  by  Pipes  et  al.  [7]  in  their  tests  with 
T300/5208  and  AS-3501-6  graphite/epoxy  laminates.  It  is 
difficult  to  discern  definite  trends  in  Rfc  with  width  ratio 
and  ply  orientation.  In  calculating  the  strength  of  loaded 
hole,  the  Rt  value  appropriate  to  the  laminate  and  hole 
configurations  should  be  used,  when  this  value  is 
unavailable,  an  approximate  value  of  Rt  must  be  used. 
Fortunately,  it  was  found  that  strength  prediction  is  not 
too  sensitive  to  the  value  of  R  .  For  example,  the  failure 
strengths  of  loaded  holes  in  T300/1Q34-C  laminates  were 
calculated  with  the  values  of  8^  *0.007,  0.018,  and  0.04  in. 
The  use  of  the  lower  and  higher  Rc  values  yielded  failure 
strengths  which  were  about  10  percent  to  20  percent 
different  from  the  one  obtained  by  the  average  Rt  value 


RAIL  SHEAR  STRENGTH,  S/S50 


VOLUME  FRACTION  OF  0°  PLIES  v0  (percent) 


Figure  19.  Variation  in  Bail  Shear  Strength  with  the 

Volume  Fraction  of  0  Degree  Plies  of  Cross  Ply 

laminates.  o  Data.  - Fit  to  Data.  S-Q 

•SOI  Voluae  Fraction  of  0  Degree  Plies  • 

19400  pai 


Figure  70.  Character i st ie  Length  in  Tension  Rt  es  Function  of  Hole  Diameter,  Width  Ratio 
and  Piy  orientaiton. 


CHARACTERISTIC  LENGTH, Rt(in) 


igore  variation  of  Char  .;rt«?r  ist  *e  length  With  Ko!@ 

Diameter.  Data  are  lor  Lata  mate  Configurations 
Qiv«n  in  Figure  20. 
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7.3)  The  Characteristic  Length  Rc 

The  value  of  Rc  was  determined  for  four  different  ply 

orientations,  as  indicated  in  Appendix  F.  The  data  are 

summarized  in  Table  4.  Each  Rc  value  in  this  table  is  the 

average  of  four  measurements. 

As  was  discussed  in  Section  6.3,  the  values  of  R  were 

c 

obtained  from  data  generated  using  loaded  holes  and  from  the 

Yamada-Sun  failure  criterion  (eq.53).  Both  the  longitudinal 

and  shear  stresses  play  a  role  in  this  criterion.  Thus, 

both  of  these  stresses  may  affect  the  value  of  R  .  The 

c 

shear  stress  has  a  significant  effect  in  those  laminates  in 

which  the  shear  stress  to  shear  strength  ratio  (o  /S)  is 

xy 

comparable  to  the  longitudinal  stress  to  longitudinal 

compressive  strength  ratio  vox/Xc).  This  situation  arises, 

for  example,  in  tQ/90].  and  £*451  laminates. 

In  calculating  Rc<  the  stresses  were  assumed  to  vary 

linearly  with  strains.  As  was  noted  previously  (Section 

6.D,  for  shear  stresses  this  assumption  cay  be  invalid* 

since  the  value  Rc  may  depend  on  the  shear  stress.  This 

assumption  say  have  affected  the  values  of  8  ,  especially 

for  the  two  cross-ply  laminates  in  Table  4.  The  effects 

introduced  in  by  the  assumption  of  linear  stress-strain 

relationship  is  unknown,*  therefore,  the  value  of  R  i»0.Q7 

c 

in.)  obtained  for  quasi -isotropic  laminates  was  adopted  in 
this  investigation. 
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Table  4  The  Characteristic  Length  in  Compression  R 

for  Fiberite  T300/1034-C.  Data  obtained  for 
D*Q.25  in,  W=2.0  in,  L=7.0  in  ,  E^l.25  in 


Ply  Orientation  Characteristic  Length  Rc  (in) 


[  (0/±45/90)3]g 

0.07 

[(902/±60/±30)2]s 

0.08 

[  ( 0/90 ) 6  3  s 

0.09 

[±45)6]s 

0.13 

SECTION  VIII 


EXPERIMENTAL  VALIDATION  OF  THE  MODEL 

In  this  section,  comparisons  between  data  and  the 
results  of  the  model  are  presented.  The  data  used  in  the 
comparisons  were  generated  during  the  course  of  this 
investigation  with  Fiberite  T300/1034-C  graphite/epoxy 
composite  having  different  geometries  and  different  ply 
orientations.  The  failure  strength  and  the  failure  modes 
were  measured  with  composites  containing  either  one  pin- 
loaded  hole  or  two  pin-loaded  holes  in  parallel,  or  two  pin- 
loaded  holes  in  series.  The  experimental  results  are 
presented  in  Figures  22  through  29. 

To  facilitate  comparisons  between  the  data  and  the 
results  of  the  model,  the  ordinates  in  these  figures 
represent  the  bearing  strength  Pg.  For  laminates  with  a 
single  hole  or  with  two  holes  in  series,  the  bearing 
strength  is  expressed  as  Pg=  P/DH.  For  laminates  containing 
two  holes  in  parallel,  the  bearing  strength  is  taken  as  PB» 
P/2DH.  P  is  the  failure  load  and  DH  represents  the  cross 
sectional  area  of  the  hole.  In  Figures  22-29  the  measured 
bearing  strengths  and  failure  modes  are  represented  by 
different  symbols. 

Tne  bearing  strengths  and  failure  modes  were  also 
calculated  by  the  model.  The  numerical  calculations  were 
performed  using  the  material  properties  listed  in  Table  3. 
The  numerical  results  are  included  in  these  figures.  The 
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calculated  bearing  strengths  are  given  by  solid  lines.  The 
calculated  failure  modes  were  not  identified  separately  as 
long  as  they  were  the  same  as  those  given  by  the  data.  In 
those  cases  where  the  calculated  failure  model  differs  from 
the  data,  the  calculated  failure  mode  is  identified  by  the 
letters  T,  B,  or  S,  next  to  the  corresponding  data  point. 
These  letters  represent  failure  in  tension,  bearing,  and 
shearout  modes. 

As  indicated  in  Figures  22-29,  for  [ (0/±45/90) and 

J  s 

[ (902/±60/±30)2]  laminates  the  calculated  failure  strengths 
agree  with  the  data  within  10  percent  to  30  percent.  The 
specimen  geometry  (hole  diameter,  edge  distance,  and  width) 
has  little  effect  on  the  accuracy  of  the  model. 

For  cross-ply  laminates( [ 0/90 3 ^  and  [±45]  )  the 
difference  between  the  calculated  bearing  strengths  and  the 
data  ranges  from  about  10  to  40  percent.  The  accuracy  is 
better  for  smaller  holes  (10  percent  for  D=  1/8  in)  and 
decreases  a  the  hole  size  increases.  The  differences 
between  the  calculated  and  measured  bearing  strengths  become 
about  40  percent  for  1/2  in  diameter  holes.  In  all  cases, 
the  calculated  values  are  conservative  and  underestimate  the 
actual  bearing  strengths.  The  reason  for  the  lower  accuracy 
of  the  model  for  cross  ply  laminates  is  most  likely  due  to 
the  assumption  that  the  shear  stress  is  linearly 
proportional  to  the  shear  strain.  Since  shear  stresses  are 
important  in  determining  the  failure  strengths  of  cross  ply 
laminates  (Section  VII),  the  use  of  nonlinear  shear  stress- 
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strain  relationships  should  improve  the  accuracy  of. the 
model  for  such  laminates. 

The  results  in  Figures  22-29  show  that  the  model  predicts 
the  failure  mode  with  good  accuracy.  Of  the  83  specimen 
configurations  tested,  the  model  failed  to  predict 
accurately  the  failure  mode  only  in  9  cases  -  -  these  cases 
being  indicated  by  the  letters,  T,  B,  or  S,  in  Figures 
22-29.  In  3  of  those  cases  where  the  model  gave  different 
failure  modes  than  the  data,  the  data  were  ambiguous. 

Failure,  in  fact,  may  have  occurred  by  the  combination  of 
two  different  modes. 

The  results  discussed  in  the  foregoing,  and  represented 
in  Figures  22-29,  show  that  the  model  provides  the  failure 
strengths  and  failure  modes  of  loaded  holes  with  reasonable 
accuracy.  The  accuracy  of  the  present  model  could  be 
improved  further  if,  instead  of  the  average  values  of  S,  Rt 
and  Rc,  the  values  corresponding  to  the  specific  geometry 
and  laminate  configuration  were  used  in  the  calculations. 

It  is  worthwhile  to  compare  the  accuracy  of  the  present 
model  with  the  accuracy  of  the  models  developed  by  previous 
investigators.  A  summary  of  the  accuracies  of  the  various 
models  is  presented  in  Table  5. 

The  accuracy  may  depend  on  the  geometry,  ply 
orientation,  and  material  properties.  Therefore,  the 
accuracies  in  Table  5  must  be  viewed  with  caution. 
Nevertheless,  the  numbers  in  this  provide  an  estimate  of  the 
magnitudes  of  error  of  the  different  models.  The  present 


model  appears  to  be  more  accurate  than  any  of  the  other 
models. 

Two  points  are  worth  noting:  First,  the  models 
developed  previously  apply  only  to  laminates  containing  a 
single  hole.  None  of  the  models  except  the  present  one 
applies  to  laminates  containing  two  holes.  Second,  of  the 
existing  models,  only  the  present  one  and  the  one  by  Garbo 
and  Ogonowski  [6]  have  been  supplemented  with  "user 
friendly"  computer  codes.  Therefore,  presently,  only  these 
two  models  can  be  used  readily.  Furthermore,  the  Garbo  and 
Ogonowski  model  yields  the  failure  strength,  but  does  not 
provide  the  mode  of  failure. 
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Figure  22, 


Bearing  Strengths  of  Fiberite  T300/J034-C 
Laminates  Containing  a  Single  Loaded  Hole. 
Comparisons  Between  the  Data  and  the  Results  of 
the  Model.  The  Failure  Modes  Calculated  by  the 
Model  are  the  Same  as  Those  of  the  Data  Unless 
Indicated  by  a  Letter  in  Parentheses  next  to 
the  Data  Point. 
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Bearing  Strengths  of  Fibcrite  T300/1034-C 
Laminates  Containing  a  Single  Loaded  Hole. 
Comparisons  Between  the  Data  and  the  Results  of 
the  Modle.  The  Failure  Modes  Calculated  by  the 
Modle  are  the  Same  as  Those  of  the  Data  Unless 
Indicated  by  a  Letter  in  Parentheses  next  to 
the  Data  Point. 
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Figure  24.  Bearing  Strengths  of  Piberite  T300/1034-C  Laminates  Containing  a  Single  Loaded 
Hole.  Comparisons  Between  the  Data  and  the  Results  of  the  Model.  The  Failure 
Modes  Calculated  by  the  Model  are  the  Same  as  Those  of  the  Data  Unless 
Indicated  by  a  Letter  in  Parentheses  next  to  the  Data  Point. 
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Figure  27.  Bearing  Strengths  of  Fiberite  T300/1034-C 

Laminates  Containing  Two  Loaded  Holes  in  Series. 
Comparisons  Between  the  data  and  the  Results  of 
the  Model.  The  Failure  Modes  Calculated  by  the 
Model  are  the  Same  as  Those  of  the  Data  Unless 
Indicated  by  a  Letter  in  Parentheses  nest  to 
the  Data  Point. 
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Figure  28.  Bearinq  Strengths  of  Fiberite  T300/1034-C 

Laminates  Containing  Two  Loaded  Holes  in  Series 
Comparisons  Between  the  Data  and  the  Results  of 
the  Model.  The  Failure  Modes  Calculated  by  the 
Model  are  the  Same  as  Those  of  the  Data  Unless 
Indicated  by  a  Letter  in  Parentheses  nest  to 
the  Data  Point . 
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Figure  79.  Searing  Strengths  of  Fibecite  T300/I034-C  Laminates  Containing  Two  Loaded 
Holes  in  Serle.  Comparisons  Between  the  Data  and  the  Results  of  the  Model 
The  Failure  Modes  Calculated  by  the  Model  are  the  Same  as  Those  of  the  Data 
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SECTION  IX 


DESIGN  CONSIDERATIONS 

As  illustrated  by  the  sample  computer  input-output  in 
Appendix  E,  the  model,  together  with  the  computer  code,  can 
readily  be  used  to  calculate  the  failure  strengths  and 
failure  modes  of  laminates  containing  a  single  pin-loaded 
hole,  two  pin-loaded  holes  in  parallel,  or  two  pin-loaded 
holes  in  series.  The  model  can  also  be  used  to  design 
joints  containing  many  pin  loaded  holes.  In  joint  design, 
it  is  desired  to  determine  the  number  of  holes,  the  hole 
diameter,  and  the  hole  positions  which  result  in  the  maximum 
failure  load  PM  and  the  maximum  failure  load  per  unit  weight 
P*M.  The  failure  load  per  unit  weight  is  defined  as 

P*  »  P/w  (59) 

where  P  is  the  failure  load,  and  w  is  the  combined  weight  of 
the  composite  w  and  the  pin  w 

w  5 

w  =  w„  +  w_  (60) 

c  s 

In  this  section,  procedures  suitable  for  calculating  PM 
and  P  M  are  illustrated  via  two  sample  problems.  In  these 
problems,  the  failure  load  of  24-ply  (thickness  H-0.125  in) 

[ ( 0/±45/90 ) 3 3 s  Fiberite  T300/1034-C  graphite-epoxy 
composites  are  determined.  The  material  properties  used  in 
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the  calculations  are  listed  in  Table  1.  The  density  of  the 
composite  is  pc  »  0.00194  lbm/'in3.  The  pin  or  pins  are 
assumed  to  be  3/4  in  long  and  to  be  made  of  steel  (density 
pg*  0.0093  lbm/in“). 

The  calculation  procedures  are  presented  in  Section  9.3 
for  joints  containing  one  or  two  holes,  and  in  Section  9.4 
for  joints  containing  three  or  more  holes.  First,  however, 
interferences  between  two  adjacent  holes,  between  the  edge 
and  an  adjacent  hole,  and  between  the  side  and  an  adjacent 
hole  are  discussed. 

9.1)  Interaction  Coefficients 


It  is  desired  to  know  under  what  conditions,  if  any,  the 
proximity  of  two  holes,  or  the  proximity  of  a  hole  to  the 
edge  or  to  the  side  of  the  laminate,  affects  the  failure 
load.  The  interaction  between  two  holes,  between  a  hole  and 
the  edge,  and  between  a  hole  and  the  side,  can  best  be 
evaluated  by  the  use  of  interaction  coefficients. 

Two  Holes  in  Parallel  The  parallel  hole  interaction 
coefficient  gH  is  defined  as 

9H  "  V(PH/2)  (61) 

Where  Pg  is  the  failure  load  of  a  GH  wide  laminate 
containing  a  single  hole,  and  PH  is  the  failure  load  of:  a 
2Gh  wide  laminate  containing  two  loaded  holes  separated  by  a 


distance  GH  (Figure  30).  When  GH  becomes  large,  the 
interaction  between  two  holes  becomes  small  (PH/2  -►  Pg)  and 
the  interaction  coefficient  approaches  unity  (gH  -►  1). 

Two  Holes  in  Series  The  series  hole  interaction 
coefficient-  gv  is  defined  as 

9V  *  PV/PT  (62) 

where  Pv  is  the  failure  load  of  a  laminate  (width  W) 
containing  two  loaded  holes  separated  by  a  distance  Gy.  PT 
is  the  failure  load  of  a  laminate  with  the  same  width 
containing  two  holes;  one  located  at  a  distance  E  from  the 
edge,  and  the  other  located  at  the  center  of  the  laminate 
(Figure  31).  When  the  hole  distance  increases,  the 
influence  of  one  hole  on  the  other  becomes  small;  the 
failure  load  Pv  approaches  PT(PV  -*•  P^)  and  gv  approaches 
unity  (gv  *  1 ) . 

Edge  Interaction  The  edge  interaction  coefficient  gE  is 
defined  as 

gE  -  VPC  (63) 

where  P<.  is  the  failure  load  of  laminates  (width  W)  with  a 
single  loaded  hole  at  distance  E  from  the  edge.  Pc  is  the 
failure  load  of  a  laminate  of  width  W  with  a  hole  in  the 
center  (Figure  32).  The  influence  of  the  edge  on  the 


failure  load  becomes  smaller  as  the  edge  distance  increases. 


When  the  hole  is  moved  to  the  center  (E-L/2),  Pg  becomes  Pc 
and  the  interaction  coefficient  becomes  unity  (gE  ♦  1). 

Side  Interaction  Coefficient  The  side  interaction 
coefficient  gg  is  defined  as 

gs  -  t>H/po  <64> 

where  P^  is  the  failure  load  of  a  laminate  (width  W) 
containing  two  loaded  holes  separated  by  a  distance  W/2.  PG 
is  the  failure  load  of  a  laminate  with  the  same-  width 
containing  two  loaded  holes  separated  by  a  distance  GH  (GH  £ 
w/2,  Figure  33).  As  the  distance  Q  between  the  side  and  the 
hole  increases  PG  *  PH  and  the  interaction  coefficient  gg 
approaches  unity. 

9*2)  Numerical  values  of  the  Interaction  Coefficients 


In  order  to  illustrate  the  trend  in  the  interaction 
coefficients  these  coefficients  were  calculated  for  Fiberite 
T300/1Q34-C  graphite-epoxy  composite  laminates  with  ply 
orientations  of  t  C 0/±45/90 ) 3 >  3 s  and  l (02/±45)3ls.  The 
results,  obtained  using  the  computer  code,  are  presented  in 
Figures  30-33.  The  most  significant  feature  of  these 
results  is  that  the  failure  load  is  not  affected 
significantly 
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a)  by  the  proximity  of  two  holes  in  parallel  when  the 
distance  between  two  holes  is  larger  than  30 

b)  by  the  proximity  of  two  holes  in  series  when  the  distance 
between  two  holes  is  larger  than  20 

c)  by  the  edge  when  the  distance  between  the  edge  and  the 
hole  is  greater  than  30. 

d)  by  the  proximity  of  side  when  the  distance  between  the 
hole  and  the  side  is  larger  than  20 

Mathematically,  these  conditions  can  be  expressed  as 


gH  -  1  PH  *  2Ps 

9y  -  *  Py  P<y 

9e  -  1  PS  pc 

9S  *  1  PG  "*  PH 


as  Gh/D  ^  3  ' 
as  Gy/D  2  2 
as  E/D  2  3 
as  Q/D  5  2 


for  l (0/145/90) L 
s 

[<02/s45)3)s 


165) 


It  is  emphasised  that  the  conditions  expressed  by  the 
above  equations  (eq.  6S )  may  not  apply  for  every  ply 
orientation.  The  conditions  at  which  the  different 
coefficients  become  unity  must  be  evaluated  separately  for 
each  laminate  lay  up. 
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Figure  10.  Interaction  Coefficient  for  Two  Holes  in  Parallel.  Results  of  the  Model. 


HOLE  DISTANCE  RATIO,Gw/D  HOLE  DISTANCE  RATIO,Gv/D 


Figure  32.  edge  Interaction  Coefficient.  Results  of  the  Model 
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Figure  33.  side  Interaction  Coefficient.  Results  of  the  Model. 
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9.3)  Laminates  with  One  or  Two  Holes 


In  this  subsection,  a  procedure  is  described  which  can 
be  used  to  size  a  laminate  containing  either  one  or  two  pin*' 
loaded  holes. 

We  consider  a  laminate  of  known  width  (w»i  in),  length 
U>Q  in)  and  thickness  (H*  0.125  in  ).  The  laminate  may 
contain  either  one  pin-loaded  hole  or  two  pin-loaded  holes 
in  parallel  or  in  series,  as  illustrated  in  Figure  34. 

It  is  desired  to  find  the  number  of  holes  (  one  or  two 
holes),  the  hole  diameter  D,  the  edge  distance  E,  and  the 
distance  between  two  holes  G,  which  result  in  the  maximum 
failure  load  PM  and  in  the  maximum  failure  load  per  unit 
weight  P*M< 

The  calculation  proceeds  along  the  following  major 
steps: 

a)  Using  the  computer  code,  the  failure  loads  of 
a  laminate  containing  a  single-loaded  hole  are 
calculated  for  different  hole  diameters  D  and 
for  different  edge  distance  ratios  E/O. 

The  failure  load  is  plotted  versus  the  edge  ratio 
E/D  (Figure  35).  The  desired  edge  ratio 
(E/O)  is  selected. 

Here,  the  edge  ratio  E/D-3  was  selected  because  the  failure 
load  reaches  a  maximum  at  the  edge  ratio  of  about  3  and 
remains  nearly  constant  at  higher  edge  ratios.  This 
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EDGE  RATIO,  E/0 


Figure  3S. 


Failure  Load  as  a  Function  of  edge  Ratio  Cor 
Laainates  Containing  a  Single  Pin-Loaded  Holes* 
Results  of  the  Model. 
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value  (E/D)  will  also  be  used  for  joints  containing  two 
loaded  holes  in  parallel  and  two  holes  in  series.  The 
reasons  for  this  choice  of  E/D  are  as  follows;  1)  For 
parallel  holes  the  interaction  between  two  holes  has  almost 
no  effect  on  the  failure  load  (gH  1  and  ?H/2  +  Pg, 
Section  9.2  ).  Hence,  when  GH/D>3  (as  is  the  case  in  the 
present  problem),  two  parallel  holes  can  be  treated  as  two 
independent  holes.  2)  For  two  holesin  series,  the 
interaction  between  the  holes  is  unimportant,  when  Gv/D>3 
(gv  +  1  and  Pv  Pc,  Section  9.2).  Two  holes  can  be 
considered  as  two  independent  holes  sharing  part  of  the 
total  load.  Hence  the  value  of  E/D=3  is  a  suitable  choice 
for  the  present  problem  when  Gv/D>3. 

b)  Using  the  computer  code,  the  failure  loads  are 
calculated  for  different  hole  diameters,  and 

for  different  hole  separations  G  for  two  holes  in 
parallel  and  for  two  holes  in  series.  The  failure 
loads  are  plotted  as  functions  of  the  hole  distance 
ratio  G/D  (Figure  36,  top).  From  these  plots, 
the  maximum  failure  load  PM  can  be  obtained. 

For  the  problem  under  consideration,  the  maximum  failure 
load  is  5000  lb.  This  load  is  achieved  by  two  0.125in  pins 
in  parallel  separated  by  a  distance  GH»0.5  in  (G^/D  ■  4). 

c)  From  the  known  values  of  the  failure  load  P, 
the  failure  load  per  unit  weight  P*  is 
calculated  using  the  expression 
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P*  =  P/[ p_WHL  +  a  (itD2/4)  (p  L  -p  H)  3  (66) 

v  5  5  C 

where  Lg  is  the  length  of  the  pin.  The 
parameter  a=1  for  a  single  hole,  a=2  for 
two  holes.  The  failure  load  per  unit  weight 
is  plotted  as  a  function  of  G/D 
(Figure  36,  bottom). 

For  the  present  problem,  the  maximum  failure  load  per  unit 
weight  P*M  is  8000  lbf/lbf  and  occurs  with  two  0.125-in 
diameter  pins  separated  by  a  horizontal  distance  GH  *  0.5  in 
(Gh/D  *  4). 

9.4)  Lvuir.utes  with  Multiple  Holes 

This  problem  is  concerned  with  laminates  containing 
several  pin-loaded  holes  spaced  evenly,  either  in  a  single 
row  or  in  two  parallel  rows,  as  illustrated  in  Figure  37. 

The  number  of  holes  in  the  laminate  with  a  single  row  of 
holes,  or  the  number  of  columns  in  the  laminates  with  two 
rows  of  holes  is 

No  "  w/gH  (67) 

It  is  desired  to  determine  the  number  of  holes  NQ,  the 
hole  size  D,  the  positions  of  the  holes  G^  and  Gv,  and  the 
edge  distance  E,  which  result  in  the  maximum  failure  load. 
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The  model  developed  in  this  investigation  can  be  applied 
only  to  laminates  containing  either  a  single  pin-loaded 
hole,  or  two  pin-loaded  holes  in  parallel  or  in  series. 
Therefore,  the  model  can  not  be  used  directly  to  calculate 
the  failure  load  of  laminate  containing  several  holes.  The 
failure  loads  of  such  laminates  can  still  be  estimated  with 
the  use  of  the  model  by  the  procedure  described  below. 

a)  The  interaction  coefficients  gE  and  gv  are  calculated  and 
plotted  in  the  manner  described  in  Section  9.1. 

b)  The  ratios  S/D  and  Gv/D  are  selected,  which  correspond  to 
the  conditions  gE  ■*  1  and  gv  1 .  In  this  investigation, 
the  values  of  both  E/D  and  Gv/D  were  selected  to  be  3 
because  both  gE  and  gv  reach  unity  at  this  ratio.  This 
E/D  ratio  is  used  for  a  single  row  of  holes.  This  is 
also  a  reasonable  choice  for  two  rows  of  holes,  because 
the  first  row  of  holes  acts  independently  of  second  rows 
of  holes,  to  a  very  large  degree. 

c)  Values  are  assumed  for  the  number  of  holes  and  the  hole 
diameter.  The  distance  GH  is  calculated  form 

Gh  «  W/Nq  (68) 
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d)  Using  the  computer  code,  the  failure  load  is  calculated 
for  a  2Gh  wide  laminate  containing  two  loaded  holes  in 
parallel  PH,  and  for  a  GH  wide  laminate  containing  two 
loaded  holes  in  series  Py.  The  interaction  parameter  gH 
is  then  calculated  for  the  geometry  under  consideration, 
according  to  the  method  given  in  Section  9.2. 

e)  The  failure  load  is  approximated  by  the  expression 


PNq-2  +  2  Pside 


(69) 


where  PN  _2  is  the  load  carried  by  the  second  (n=2)  through 
o 

the  next  to  last  (n=NQ-1)  pins,  and  Psi(Je  is  the  load 
carried  by  the  first  (n=>1)  and  last  (n»NQ)  pins.  Thus,  the 
failure  load  of  a  laminate  containing  one  or  two  rows  of 
holes  is 


Pr1  "  < (NQ-2)/2)gHPH  +9SPH  (one  row)  (70) 

Pr2  *  (No”2^H^PV  *  2gH9SPV  (two  rows*  (71) 

If  Q*Gh/2  then  gg  is  equal  to  unity.  This  is  the  case 
in  the  present  problem.  Accordingly, 

Pr,  -  <(N0-2)/2)gHPH  *  PH  (72) 

Pr2  ‘  +  29hPV 


(73) 
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f)  The  failure  load  per  unit  weight  P*  is  calculated 

p*rl,r2  ■  Pri,r2/tl>cWKLtaNo(,,D2'/4)<‘!sLs-',cH)1  (74) 

where  the  subscripts  rl  and  r2  refer  to  one  row  or  two  rows 
of  holes,  respectively. 

g)  The  calculations  are  repeated  tor  different  values  of  NQ 
and  D.  The  failure  load  Pf1  r2  and  the  failure  load  per 
unit  weight  P*f1  f2  are  plotted  as  functions  of  Nq.  From 
these  figures,  the  maximum  failure  load  PM  and  the 
maximum  failure  load  per  unit  weight  P*M  are  determined. 

In  the  present  problem  a  4  in  wide  and  10  in  long 
composite  laminate  was  considered.  The  ply  orientation  is 
[(0/t45/±90)3]s.  The  material  properties  are  listed  in 
Table  3.  The  procedures  gave  the  maximum  failure  load  PM  • 
23400  lbf  when  there  are  twelve  0.125  in  diameter  holes 
arranged  in  two  rows  of  holes  (Figure  38).  The  maximum 
failure  load  per  unit  weight  (P*M  ■  67000  lbf /lbf )  is 
achieved  with  twelve  0.125  in  diameter  holes  in  a  single  row 
(Figure  38). 

9.5)  Failure  Mode 

The  results  generated  by  the  computer  code  also  show  the 
modes  of  the  failure.  The  changes  in  the  modes  of  failure 
with  the  number  of  holes  NQ  are  illustrated  in  Figure  39. 

In  the  present  sample  problem,  at  the  condition  of  the 
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maximum  failure  load  (NQ>12)  the  failure  mode  is  in  tension. 
Failure  in  such  mode  often  happens  quite  suddenly.  In  some 
situations  if  might  be  preferable  to  choose  a  design  in 
which  failure  occurs  by  a  less  sudden  failure  mode.  For 
example,  failure  would  have  occurred  in  bearing  mode  if,  in 
the  present  problem,  a  hole  diameter  of  0.125  in  were 
chosen,  and  the  number  of  holes  were  taken  to  be  NQ«6. 
However, this  would  have  resulted  in  a  30  percent  to  40 
percent  reduction  in  the  failure  load. 


FAILURE  LOAD  PER  UNIT  WEIGHT  FAILURE  LOAD.P  KlO's(lbf) 


NUMBER  OF  HOLES.  Nq  number  of  HOLES.  Nq 

Figure  38.  Failure  Load  (Top)  and  Failure  Load  Per  Unit 

Weight  (Bottom)  of  Laminates  Containing  One  Roe 
(Left)  and  Two  Roes  (Right)  of  Pin-Loaded  Kolas* 
Results  of  the  Model. 
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Figure  39.  Failure  Modes  of  Single  Row  (Left)  and  Two  Rows 
(Right)  of  Pin-Loaded  Holes.  Results  of  the  Model 


SECTION  2 


SUMMARY  AND  CONCLUSIONS 


The  following  major  tasks  were  completed  during  the 

course  of  this  investigation: 

a)  A  model  and  a  computer  code  were  developed  which  can  be 
used  in  the  design  of  mechanically-fastened  composite 
joints  involving  fiber  reinforced  laminates.  The  model 
can  be  used  to  determine  the  failure  loads  and  failure 
modes  of  laminates  containing  a  single  pin-loaded  hole, 
two  pin-loaded  holes  in  parallel,  and  two  pin-loaded 
holes  in  series. 

b)  Experimental  procedures  were  developed  to  determine  the 
characteristic  lengths. 

c)  Tests  were  performed  to  determine  the  values  of  the  rail- 
shear  strength  and  the  characteristic  lengths  of  Fiberite 
T300/1034-C  composites,  and  to  evaluate  the  effects  of 
geometry  and  laminate  lay  up  on  these  parameters. 

d)  A  series  of  tests  was  performed  measuring  the  failure 
strengths  and  failure  modes  of  Fiberite  T300/1034-C 
laminates  containing  a  single-pin  loaded  hole,  two  pin- 
loaded  holes  in  parallel,  and  two  pin-loaded  holes  in 
series. 


101 


102 


e)  Comparisons  were  made  between  the  data  and  the  results  of 
the  model.  Good  agreements  were  found  between  the 
analytical  ana  the  experimental  results. 

f)  Procedures  were  developed  for  the  design  of  composite 
laminates  containing  one,  two,  or  more  pin-loaded  holes. 

The  model  was  developed  on  the  basis  of  the  following 
assumptions:  a)  classical,  two-dimensional  laminate  plate 
theory,  and  b)  linear  relationship  between  the  stresses  and 
strains.  Good  agreements  between  the  results  of  the  model 
and  the  data  suggest  that  these  assumptions  are  reasonable 
for  a  wide  range  of  problems.  Three  dimensional  stress 
distributions  and  nonlinear  stress-strain  relationships 
could  oe  incorporated  into  the  model  in  the  future. 
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appendix  a 


The  components  of  the  matrix  QP.j  appearing  in  Eg.  (12)  are 


Q  ii“Q  1 ,cos  h+2(QP12+2QP66)sin2n+QP22sin4n ' 

®  1 2“ (qE>i  1  ■+‘qP22~4qP66 ^ sin2-ncc>s2Tr^Qp  12(sin4-n+cos4n) 

QP22«QPi  lSin4n+2(QP  2+2QP66)sin211cos  VqP22cos4t1 

QP13*  <qP1  l"QP12“2QP33)sinT|COs3'n  +  (QP1  2“QP22+2qP33  )sin3,ncos^ 

SP23"(QPH"Qp12'2QP33)sin3TlcoST1+(QP12-QP22+2QP33)sinncos3Ti 

^P33"  (qP1  1  ^qE>22~2qP1  2_2qP33  > s i n 2ticos2ti+QP3 3  ( si n4Ti+cos4n ) 
in  which 

0Pu  •  sVc-oVV 

Qp,2  •  <»pla*F,>/<»V1J»,ll>'  UP2,EP,/U-UPI2UP2,) 

'f22  *  EP2/(1-»P12,.P2i) 

0*3;  *  “P,a 
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The  superscript  p  denotes  the  material  properties  of  the 


p-th  ply,  and  the  angle  7}  is  measured  from  the  x«-axis  to 


*  P  P  P 

the  x-axis.  E  y  E  ^  an&  G  12  are  longitudinal , 


transverse,  and  shear  moduli  of  the  p-th  ply,  respectively. 


p  p 

V  12  an^  V  21  are  P°isson's  ratios  for  the  p-th  ply  and 


satisfy  the  relation 


P  /trp 
v  12/E  t 


P  /WP 


•  .W'  *  .* 
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APPEND I Z  B 


The  Coordinate  Transformation  matrix  T. 


i] 


[T]  « 


2 

CCS  T} 

sin2n 

simricosn 


sin2n 
COS  T) 

sinncosn 


“■ 

2sinncosn 

-2sinncos,n 

2  •  2 
cos  Tpsin  n 


The  angle  ti  is  measured  from  the  x^-axis  to  the  x  axis. 
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APPENDIX  C 

The  Finite  Element  Mesh  Generator 

The  mesh  generator  generates  306  quadrilateral  elements 
for  a  single  hole,  and  612  and  655  elements  for  two  holes  in 
parallel  and  two  holes  in  series,  respectively.  To  control 
costs,  this  number  was  held  fixed  (the  reader  should  note 
that  612  elements  involve  a  matrix  of  size  1400  x  300).  The 
mesh  is  designed  in  such  way  that  the  characteristic  curve 
can  be  encompassed  by  a  square  of  size  2z  x  2z,  in  which  a 
fine  mesh  was  generated  and  outside  of  which  the  calculated 
stresses  are  still  reasonably  accurate  (see  Figures  5,  6, 

7).  A  suitable  value  of  z  for  a  given  geometry  has  to  be 
determined  before  calculating  the  stresses.  Mathematically, 
this  problem  may  be  stated  : 

Optimize  z 

subject  to  (  Rc  ♦  D/2  )  s  z  s  W/2 

z  S  E 

Equation  (C.l)  can  be  rewritten  as 

0  S  (  z  -  R  -  D/2  )  £  (  W/2  -  R  -  D/2  )  (C.3) 

v  C 

Assume  that  E  is  large  enough  that  eq.(C.2)  is  always 
satisfied.  Assume  that  z/D  is  a  function  of  W/d  and  R„/D, 

W 


(c.l) 

(C.2) 


Ill 


I 


and  can  be  expressed  by 

z/D  »  f(  W/2D  -  Rc/D  -1/2  )  -  f(  C/D  )  (C.4) 

where  c/D  ■  W/2D  -  Rc/D  -  1/2  and  i  is  some  unknown 
function. 

Assume  that  z/D  is  a  second  order  polynomial  function  of 
C/D.  Then  eq.(C.4)  can  be  written  as 

z/D  ■  a(c/D)2  +  b( C/D)  +  c  (C.5) 

This  equation  reflects  the  general  trend  that  as  W,  and 

hence  C»  increases,  z  must  be  made  bigger,  since  the  total 

number  of  ilements  is  constant.  Three  conditions  are 

necessary  to  determine  the  constants,  a,  b,  and  c. 

When  W/2  ■  R_  ♦  D/2  (  C  •  0  ) ,  then  the  only  choice  of  z 
c 

is 

z/D  •  Rc/D  ♦  1/2  -  W/2  (C.6) 

Substituting  eq.(C.6)  into  (C.5),  gives 
v  2300 

c  •  Rc/D  ♦  1/2  (C.7) 

When  w  changes  from  w,  to  w2,  say,  (  changes  by  the 
amount  AC  ■  c,  -  C2 .  The  writer  has  found  from 
computational  experience  that  the  following  changes  in  z 


112 


-E 


generated  results  close  to  known  analytical  solutions. 

1  S  5  S  1.5  A(z/D)  e  (1/4)  A(C/D)  ( C . 8 > 

5*2  A ( z/D)  -  4  A (5/D)  (C.9) 

Imposing  these  conditions,  we  have 

at  5  »  1,  A(z/D)/A(5/D)  =d(z/D)/d(5/D)  •  2a  +  b  »  1/4 

(C. 10) 

at  5  .5,  A<z/D)/AU/D)  =d(z/D) /d(  5/D)  =  5a  +b  =  4  (C.11) 

From  eqs  (C.10)  and  (C.11),  give 

a  -  0.05  (C. 12) 

b  ■  0.15  (C.13) 

As  a  result,  eq.(C.5)  becomes 

z  m  D  (  0.05  (W/2D  -  R  /D  -  1/2)2  ♦  0.15  (W/2D  -R  /D  -  1/2) 

C  * 

♦  (R  /D  ♦  1/2))  (C.14) 

G 

Using  this  result,  excellent  agreement  between  the 
computational  results,  and  Timoshenko's  and  De  Jong's 
solutions  were  obtained  (Figures  13  and  14,  Section  v). 
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Appendix  D 


Shape  Function  Used  in  the  Finite  Element  Code 


In  the  isoparametric  element,  the  geometry  and  the 
displacement  of  the  element  are  described  in  terms  of  the 
shape  function  Na,  by  a  transformation  from  a  master  element 
in  the  r-s  coordinate  system  to  the  element  in  the  x1~x2 
coordinate  system  (Figure  40). 


X.  -  Na(Y,s)  xia  i-1.2 

ui  *  Vy,s>  ^ia  o-1, 2, 3, or  4 

N  (r # s ) - 1/4 ( 1*rr  )  (  Hss  )  -ISr.sSI 
a  a  a 


Here  xt  is  the  coordinate  of  node  a  in  the  indirection,  q? 
is  the  displacement  of  node  a  in  the  i-direction,  and  rft  and 
s  are  the  coordinates  of  node  a  referred  to  the  master 

ot 

element.  Note  the  property 


WV  “ 


1,  if  a-8 

< 

0,  if  a*B 


r 
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Element  in  x,-x2  coordinates 


Master  element  in  local 
r-s  coordinates 


Figure  40  Geometry  of  an  Element  Used  in  the  Finite  Element 
Calculations;  Left:  Element  in  the  x.-x? 
Coordinate  System.  Right:  Element  (Master 
Element)  in  the  Local  (r-s)  Coordinate  System, 
x^  is  the  Coordinate  of  Node  a  in  the  i 
Direct  ion.  qi  is  the  Displacement  of  Node 
at  in  the  i  Direction  and  (r  ,s  )  are  the 
Coordinates  of  Node  a  in  th8  r^s  Coordinate 
System,  i»l,2,  a-1,2,3,  or  4. 
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Listing  of  a  Sample  of  Input-Output  of  the  Computer  Code 
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********************************************************* 


—  <<  BOLTED  JOINTS  »  - 


THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  PREDICT 
THE  FAILURE  LOAD  AND  THE  FAILURE  MODE  OF 
BOLTED  COMPOSITE  JOINTS. 


FU-KUO  CHANG,  RICHARD  A.  SCOTT,  GEORGE  S.  SPRINGER 
MECHANICAL  ENGINEERING  AND  APPLIED  MECHANICS 
THE  UNIVERSITY  OF  MICHIGAN 
ANN  ARBOR,  MI  48109 
APRIL  30,  1983 


CAPABILITIES: 


THIS  PROGRAM  HAS  THE  CAPABILITY  TO  DEAL  WITH  THREE 
TYPES  OF  BOLTED  COMPOSITE  JOINTS  DEFINED  AS  FOLLOWS: 


TYPE  1 
TYPE  2 
TYPE  3 


JOINTS  WITH  A  SINGLE  HOLE 

JOINTS  WITH  TWO  IDENTICAL  HOLES  IN  A  ROW 

JOINTS  WITH  TWO  IDENTICAL  HOLES  IN  TANDEM 


SEE  FIGURE  BELOW : 


TYPE  1 


YPE  2 


TYPE  3 
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THIS  PROGRAM  CAN  ALSO  HANDLE  THE  FOLLOWING 
LOADING  CONDITIONS: 

(A).  PIN  OR  PINS  CARRY  ALL  THE  APPLIED  LOAD. 

IB).  PIN  OR  PINS  CARRY  ONLY  A  FRACTION  OF  THE  TOTAL 
LOAD  APPLIED  AT  THE  BOTTOM  OF  THE  JOINT.  THE 
REST  OF  THE  LOAD  IS  CARRIED  BY  THE  UPPER  END. 

SEE  FIGURE  BELOW: 


P2 


/ 


\ 


/ 


\P1 


P-P1*P2,  P>P1  AND  P>?2  OR  P-P1 

P:  THE  APPLIED  LOAD 
P 1 : LOAD  CARRIED  BY  THE  PIN  (PINS) 
P2: BY-PASSED  LOAD 


\ 


/ 


P 


FOR  EACH  TYPE  OF  JOINT,  THIS  PROGRAM  CAN  HANDLE 
THE  FOLLOWING  SITUATIONS: 


(A).  DIFFERENT  PLY  ORIENTATIONS 

(B>.  DIFFERENT  MATERIAL  PROPERTIES  (SYMMETRIC  LAMINATE) 
(C).  DIFFERENT  GEOMETRICAL  CONFIGURATIONS  INCLUDING 
DIFFERENT  HOLE  SITES,  HOLE  POSITIONS,  JOINT 
THICKNESSES,  AND  JOINT  LENGTHS , 

.  RESTRICTIONS: 


THE  PROGRAM  IS  BASED  ON  THE  FOLLOWING  ASSUMPTIONS: 

(1) .  A  UNIFORM  TENSILE  LOAD  IS  APPLIED  SYMMETRICALLY 

WITH  RESPECT  TO  THE  CENTS KLINE  OF  THE  PLATE. 

(2) .  THE  LAMINATE  IS  SYMMETRIC. 
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(3) .  HOLE  SIZES  ARE  EQUAL  IN  EACH  JOINT  WITH  TWO  HOLE 

(4) .  PIN  IS  RIGID.  THE  PIN  SUPPORT  IS  ALSO  RIGID. 


ANALYSIS: 


THE  STRESSES  ARE  CALCULATED  USING  A  FINITE  ELEMENT  METHOD 
FORMULATED  ON  THE  BASIS  OF  TWO  DIMENSIONAL  CLASSICAL 
LAMINATION  PLATE  THEORY.  THE  FAILURE  LOAD  AND  FAILURE 
MODE  ARE  CALCULATED  USING  THE  CHANG-SCOTT-SPRI NGER  FAILURE 
HYPOTHESIS  TOGETHER  WITH  THE  YAMADA-SUN  FAILURE  CRITERION 


INPUT  INSTRUCTIONS 


ENTER  MATERIAL  PROPERTIES  ***• 


DO  YOU  WANT  TO  USE  GRAPHITE/EPOXY  T30Q/1Q34-C? 

ENTER  YES  OR  NO 

yes 

MATERIAL  PROPERTIES  OF  T300/1034-C  : 

PSI 
PS  I 
PSI 

PSI 
PSI 
PSI 

CHARACTERISTIC  LENGTH  (TENSION):  0.0180  INCH 

CHARACTERISTIC  LENGTH (COMPRESS I ON) :  0.0700  INCH 


LONGITUDINAL  YOUNGS  MODULUS: 
TRANSVERSE  YOUNGS  MODULUS: 

SHEAR  MODULUS: 

POISSON  RATIO: 

LONGITUDINAL  TENSILE  STRENGTH: 
LONGITUDINAL  COMPRESSIVE  STRENGTH: 
LAMINATE  SHEAR  STRENGTH: 


21300000.00000 
1700000.00000 
397000.00000 
0.30000 
251000. 00Q00 
200000.00000 
1S400. 00000 


JOINT  TYPE  SELECTION 

TYPE  1  :  JOINT  WITH  A  SINGLE  HOLE 
TYPE  2  :  JOINT  WITH  TWO  HOLES  IN  ROW 
TYPE  3  :  JOINT  WITH  TWO  HOLES  IN  TANDEM 

WHICH  TYPE  OF  JOINT  DO  YOU  WANT  TO  SELECT? 
ENTER  1,  2.  OR  3. 

1 


1  1  9 


DO  YOU  CONSIDER  A  BY-PASS  LOAD? 

ENTER  YES  OR  NO 

no 

THE  FOLLOWING  GEOMETRIC  PARAMTERS  MUST  BE  SPECIFIED: 

(A)  DIAMETER  OF  THE  HOLE,  D 

<D  SHOULD  BE  LESS  THAN  1  INCH  FOR  DEPENDABLE  RESULTS) 

(B)  WIDTH  OF  THE  JOINT,  W 

(C)  LENGTH  OF  THE  JOINT,  L 

(D)  EDGE  DISTANCE  OF  THE  JOINT,  E 

(E)  DISTANCE  BETWEEN  THE  CENTERS  OF 
TWO  HOLES,  S 


SEE  FIGURE  BELOW: 


I-  S  -|  |  | 

*  v  — 

S 


THE  DIAMETER  MUST  BE  IMPUTED  IN  INCHES,  THE  OTHER 
GEOMETRIC  PARAMETERS  MAY  BE  EITHER  IN  INCHES 
OR  AS  A  RATIO  TO  DIAMETER  (PARAMETER/  DIAMETER) 


ENTER  THE  HOLE  DIAMETER  IN  INCHES 
0.25 

DC  YOU  WISH  TO  ENTER  ALL  GEOMETRIC  PARAMETERS 
IN  TERMS  OF  DIAMETER  RATIO  (PARAMETER/DIAMETER)  ? 
ENTER  YES  OR  NO 
Y 


ENTER  THE  WIDTH  TO  DIAMETER  RATiO: 


8 
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ENTER  THE  EDGE  TO  DIAMETER  RATIO: 

3 

ENTER  THE  LENGTH  TO  DIAMETER  RATIO: 


20 


INPUT  THE  JOINT  THICKNESS  AND  THE  PLY  ORIENTATIONS 
(SYMMTRIC  LAMINATE  ONLY) 

THE  PLY  ORIENTATION  AND  THE  NUMBERS  OF  PLIES  IN  THE 
PLY  GROUP  HAVE  TO  BE  SPECIFIED. 

*1,  THE  PLY  GROUP  IS  DEFINED  AS  A  GROUP  OF  PLIES 
HAVING  THE  SAME  PLY  ORIENTATION. 

*2.  EACH  PLY  ORIENTATION  IS  MEASURED  FROM  THE 
LOADING  DIRECTION  TO  THE  FIBER  DIRECTION. 

THE  ANGLE  IS  POSITIVE  CLOCKWISE  AND 
NEGATIVE  COUNTERCLOCKWISE 

SEE  FIGURE  BELOW 


ANGLE/ 

♦  / 

i 

/ 

/ 

/ 

/ 


*  HOLE 


ENTER  THE  JOINT  THICKNESS  IN  INCHES 
0.125 

ENTER  THE  TOTAL  NUMBER  OF  PLY  GROUPS ♦ 

•••  THE  MAXIMUM  NUMBER  OF  THE  PLY  GROUPS  <  100 
INPUT  AN  INTEGER 


4 

ENTER  THE  PLY  ORIENTATION  OF  EACH  PLY  GROUP 
ENTER  THE  PLY  ORIENTATION  OF  PLY  GROUP  ? 


1 21 


IN  DEGREES 
0 

ENTER  THE  NUMBER  OF  PLIES  IN  PLY  GROUP  t  IN  INTEGER 
6 

ENTER  THE  PLY  ORIENTATION  OF  PLY  GROUP  2 
IN  DEGREES 

45 

ENTER  THE  NUMBER  OF  PLIES  IN  PLY  GROUP  2  IN  INTEGER 
6 

ENTER  THE  PLY  ORIENTATION  OF  PLY  GROUP  3 
IN  DEGREES 

-45 

ENTER  THE  NUMBER  OF  PLIES  IN  PLY  GROUP  3  IN  INTEGER 
6 

ENTER  THE  PLY  ORIENTATION  OF  PLY  CROC P  4 

IN  DEGREES 


90 

ENTER  THE  NUMBER  OF  PLIES  IN  PLY  GROUP  4  IN  INTEGER 
6 

DO  YOU  WANT  TO  HAVE  A  LIST  OF  THE  INPUT  DATA  ? 

ENTER  YES  OR  NO 

Y 

LIST  OF  DATA 
JOINT  TYPE  SELECT! ON*  ! 

LOAD  TYPE  SELECTION:  0.0  *  0?  BY-PASSED  LOAD 

*  GEOMETRY  »:  (INCHES) 

DIAMETER  WIDTH  EDGE  THICKNESS  LENGTH 

0.7500  0.1250  5.0000 


0.2500 


2. 0000 
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<  GROUP  ORIENTATION  >: 


TOTAL 

PLY 

GROUP  NO. *  4 

GROUP 

1 

ORIENTATION®  0.0 

THICKNESS® 

0.03125 

INCH 

GROUP 

2 

ORIENTATION®  45.000 

THICKNESS® 

0.03125 

INCH 

GROUP 

3 

ORIENTATIQN®-45. 000 

THICKNESS® 

0.03125 

INCH 

GROUP 

4 

ORIENTATION*®  90.000 

THICKNESS® 

0.03125 

INCH 

<  MATERIAL  PROPERTIES  >  : 


LONGITUDINAL  YOUNGS  MODULUS  ; 
TRANSVERSE  YOUNGS  MODULUS  : 

SHEAR  MODULUS: 

POISSON  RATIO: 

LONGITUDINAL  TENSILE  STRENGTH: 
LONGITUDINAL  COMPRESSIVE  STRENGTH; 
LAMINATE  SHEAR  STRENGTH: 


21300000. 00000  PSI 
1700000.00000  PSI 
897000.00000  PSI 
0.30000 

251000.00000  PSI 
200000.00000  PSI 
19400.00000  PSI 


CHARACTERISTIC  LENGTH  (TENSION):  0.0180  INCH 

CHARACTERISTIC  LENGTH (COMPRESS I ON ) ;  0.0700  INCH 


DO  YOU  WANT  TO  MAKE  ANY  CHANGE  IN  YOUR  DATA? 
ENTER  YES  OR  NO 


no 


*  THB  STRENGTH  PREDICTION  OF  FASTENED  COMPOSITE  JOINTS  > 


LIST  OF 

INPUT 

JOINT  TYPE 

SELECTION* 

1 

LOAD  TYPE  S 

■ELECTION: 

0,0  «  OF 

SS “PASSED 

LOAD 

GEOMETRY 

1  >  i  (INCHES) 

OI A METER 

WIDTH 

EDGE 

THICKNESS 

LENGTH 

0.2 S0Q 

2.0000 

0.7500 

0, 1250 

S . 0000 

"W't 


7-CC-'  > ’ :  D  jA  '♦vH  - 

-  -A-"-  uk.' 
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:  GROUP  ORIENTATION  >  : 

TOTAL  PLY  GROUP  NO.-  4 

GROUP  1  ORIENTA?ION=  0.0  THICKNESS-  0.03125  INCH 

GROUP  2  ORIENTATION®  45.000  THICKNESS-  0.03125  INCH 

GROUP  3  ORIENTATION- -4 5. 000  THICKNESS-  0.03125  INCH 

GROUP  4  ORIENTATION-  90.000  THICKNESS-  0.03125  INCH 


MATERIAL  PROPERTIES: 

LONGITUDINAL  YOUNGS  MODULUS: 
TRANSVERSE  YOUNGS  MODULUS: 

SHEAR  MODI ’.US: 

POISSON  RATIO: 

LONGITUDINAL  TENSILE  STRENGTH: 
LONGITUDINAL  COMPRESSIVE  STRENGTH: 
LAMINATE  SHEAR  STRENGTH: 


21300000.00000  PSI 
1700000.00000  PSI 
897000.00000  PSI 
0.30000 

251000.00000  PSI 
200000.00000  PSI 
19400.00000  PSI 


CHARACTERISTIC  LENGTH  (TENSION):  0.0180  INCH 

CHARACTERISTIC  LENGTH (COMPRESS I ON ) :  0.0700  INCH 


LIST  OF  OUTPUT 


«  FAILURE  LOAD  AND  FAILURE  MODE  > 


THE  MAXIMUM  LOAD  (  P  )• 

THE  BEARING  STRENGTH (  P/(D«H)  )• 
(N  :  THE  LAMINATE  THICKNESS) 


3012.7  LB 
90406.5  PSI 


THE  FAILURE  MODE  :  BEARING  MODE,  AT  THE  ANGLE  8.43?  DEGREE 


THE  FAILURE  ANGLE  IS  DEFINED  IN  THE  FOLLOWING  FIGURE  : 


t:  *  ^ - . ~i  ■‘-—i  T . .  '  " 


‘F’  -  v 

J  *>  -  > 


*  / 

ANGLE 

/ 

\ 

/ 

1 

\  ♦ 

*  / 

\ 

/ 

*  / 

*  * 

/ 

HOLE 

1 

* 

HOLE 

*  THE  INITIAL  FAILED  PLY  GROUP  (AT  THE  MAXIMUM  LOAD)  • 
THE  PLY  ORIENTATION  OF  THIS  PLY  GROUP*  0.0 


DO  YOU  WANT  TO  RUN  THE  PROGRAM  AGAIN? 
ENTER  YES  OR  NO 

y 


DO  YOU  WANT  TO  MAKE  ANY  CHANGE  IN  YOUR  DATA? 
ENTER  YES  OR  NO 


y 

WHICH  PART  OF  THE  DATA  DO  YOU  WANT  TO  CHANGE  ? 

( 1 )  .  JOINT  TYPE  AND  GEOMETRY . 

(2) .  PLY  ORIENTATION. 

(3) .  MATERIAL  PROPERTIES. 

ENTER  1  *  2  <  OR  3. 

1 

JOINT  TYPE  SELECTION 


TYPE  I  :  JOINT  WITH  A  SINGLE  HOLE 
TYPE  2  :  JOINT  WITH  ".WO  HOLES  IN  ROW 
TYPE  3  :  JOINT  WITH  TWO  HOLES  IN  TANDEM 


WHICH  TYPE  OF  JOINT  DO  YOU  WANT  TO  SELECT? 
ENTER  1.  2.  OR  3. 
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3 


DO  YOU  CONSIDER  A  BY-PASS  LOAD? 

ENTER  YES  OR  NO 

n 

THE  FOLLOWING  GEOMETRIC  PARAMTERS  MUST  BE  SPECIFIED: 

(A)  DIAMETER  OF  THE  HOLE,  D 

(D  SHOULD  BE  LESS  THAN  1  INCH  FOR  DEPENDABLE  RESULTS) 

(B)  WIDTH  OF  THE  JOINT,  W 

(C)  LENGTH  OF  THE  JOINT,  L 

(D)  EDGE  DISTANCE  OF  THE  JOINT,  E 

(E)  DISTANCE  BETWEEN  THE  CENTERS  OF 
TWO  HOLES,  S 


SEE  FIGURE  BELOW: 


E 

I’  S  -|  |  | 

•  * 

S 

L  • 


THE  DIAMETER  MUST  BE  INPUTED  IN  INCHES,  THE  OTHER 
GEOMETRIC  PARAMETERS  MAY  BE  EITHER  IN  INCHES 
OR  AS  A  RATIO  TO  DIAMETER  (PARAMETER/  DIAMETER) 


ENTER  THE  HOLE  DIAR:  .ES  IN  INCHES 


0.25 

DO  YOU  WISH  TO  ENTER  ALL  GEOMETRIC  PARAMETERS 
IN  TERMS  OP  DIAMETER  RATIO  i PARAMETER/DI AMETER )  ? 
ENTER  YET  OR  NO 


y 


ENTER  THE  WIDTH  TO  DIAMETER  RATIO 


8 

ENTER  THE  EDGE  TO  DIAMETER  RATIO: 

3 

ENTER  THE  LENGTH  TO  DIAMETER  RATIO: 

20 

ENTER  THE  TWO  HOLE  DISTANCE  TO  DIAMETER  RATIO: 
3 

DO  YOU  WANT  TO  HAVE  A  LIST  OF  THE  INPUT  DATA  ? 
ENTER  YES  OR  NO 

n 


<  THE  STRENGTH  PREDICTION  OF  FASTENED  COMPOSITE  .301  NTS 


LIST  OF  INPUT 


JOINT  TYPE  SELECTION-  3 

LOAD  TYPE  SELECTION:  0.0  X  OF  BY-PASSED  LOAD 

«  GEOMETRY  >  :  (INCHES) 

DIAMETER  WIDTH  EDGE  THICKNESS  LENGTH 

0.2500  2.0000  0.7500  0.1250  5.0000 

DISTANCE  BETWEEN  THE  TWO  HOLES  (INCHES) 

0.7500 


<  GROUP  ORIENTATION  >  : 
TOTAL  PLY  GROUP  NO.-  4 

GROUP  1  ORIENTATION-  0.0 

GROUP  2  ORIENTATION-  45.000 

GROUP  3  ORIENTATION— 45.000 


THICKNESS- 

THICKNESS- 

THICKNESS" 


0.03125  INCH 
0.Q312S  INCH 
0.03125  INCH 


127 


fc 

£ 


GROUP  4  ORIENTATION-  90.000  THICKNESS*5  0.03125  INCH 


MATERIAL  PROPERTIES: 


LONGITUDINAL  YOUNGS  MODULUS: 
TRANSVERSE  YOUNGS  MODULUS: 

SHEAR  MODULUS: 

POISSON  RATIO: 

LONGITUDINAL  TENSILE  STRENGTH: 
LONGITUDINAL  COMPRESSIVE  STRENGTH: 
LAMINATE  SHEAR  STRENGTH: 


21300000.00000 

1700000.00000 

897000.00000 

0.30000 

251000.00000 

200000.00000 

19400.00000 


PSI 

PSI 

PSI 

PSI 

PSI 

PSI 


CHARACTERISTIC  LENGTH  (TENSION):  0.0180  INCH 

CHARACTERISTIC  LENGTH (COMPRESS I ON ) :  0.0700  INCH 


LIST  OF  OUTPUT 


<  FAILURE  LOAD  AND  FAILURE  MODE  > 


THE 

MAXIMUM 

LOAD  (  P  )• 

5346.7 

LB 

THE 

BEARING 

STRENGTH (  P/(D*H)  )• 

171093.2 

PSI 

(H  :  THE  LAMINATE  THICKNESS) 


THE  FAILURE  MODE  •  SHEAROUT  MODE.  AT  THE  ANGLE  47.812 
DEGREE 


THE  FAILURE  ANGLE  IS  DEFINED  IN  THE  FOLLOWING  FIGURE  : 


•  / 

ANGLE 

/ 

\ 

/ 

• 

\  * 

•  / 

\ 

/ 

♦  / 

• 

/ 

HOLE 

* 

HOLE 

*  THE  INITIAL  FAILED  PLY  GROUP  (AT  THE  MAXIMUM  LOAD) 
THE  PLY  ORIENTATION  OF  THIS  PLY  GROUP*  45.000 

*««  THE  FAILURE  INITIATED  FROM  THE  BOTTOM  HOLE 


LOAD  CARRIED  BY  THE  TOP  PIN  *  2028.644466  LB 

LOAD  CARRIED  BY  THE  BOTTOM  PIN  *  3318.018795  LB 


DO  YOU  WANT  TO  RUN  THE  PROGRAM  AGAIN? 
ENTER  YES  OR  NO 

no 

#Execution  terminated 
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APPENDIX  P 


Summary  of  Data  for  Calculating  R.  and  k 

w  C 

This  Appendix  contains  the  data  which  were  generated  to 
determine  the  rail  shear  strength  S  and  the  characteristic 
lengths  and  Rc  for  Fiberite  T300/1Q34-C  graphite  epoxy 
laminates. 

Notations  used  in  Tables  6-14 


D  hole  diameter 

w  specimen  width 

L  specimen  length 

H  specimen  thickness 

P  edge  distance 

P  failure  load  under  tension 

pavg  »v«rage  failure  load  under  tension 

S  rail  shear  strength 

Savg  average  rail  shear  strength 

R.  characteristic  length  in  tension 

&c  characteristic  length  in  compression 


( in) 
(in) 

( in ) 

( in) 
(in) 
(lbf) 
(lbf ) 
(psi ) 
(psi) 
( in) 

( in) 


130 


O' 

m 

03 

V0 

> 

00 

(3 

© 

(0 

* 

03 

(N 

«1 

O' 

CM 

n 

H 

0“. 

0-4 

« 

i  « 

* 

VC  O 

o 

o 

o 

o  o  o 

O  -s  O 

0 

O'  r» 

(N 

m 

o 

o  in  * 

«*  o»  in 

44 

W1 

O  00 

C*> 

(V) 

in 

in  m  o» 

r~  in  m 

0 

O'  O' 

CO 

m 

CO 

NON 

^  n« 

» 

c 

iH 

0-4 

M 

rH 

r-4  r-4  H 

H  HH 

0 

O' 

e 

0 

(0 

in 

cn 

m 

u 

a 

CM 

o 

o 

0 

X 

0-4 

r-4 

H 

> 

UI 

t 

4 

** 

0 

o 

O 

o 

o 

V 

ov 

X 

\ 

in  in 

O 

o 

o 

o 

o 

44 

o 

r*' 

o 

o 

o 

o 

o 

W 

a 

•  • 

• 

• 

• 

• 

* 

c 

4^ 

r-  r* 

CD 

CD 

03 

CO 

CD 

*•4 

>4  in 

*~i  1) 

! 

■o 

CL  X 

lT> 

O 

O 

O  lf»L3 

o  in  in 

0 

u 

W  lf> 

f- 

o 

o 

o  i>  r~ 

o  r*  r- 

■o 

<n  c 

X 

•  * 

• 

• 

• 

i  i  « 

i»» 

3 

w  •*« 

H  H 

m 4 

cm 

CM 

f'i  H  H 

(N  H  H 

H 

o 

u 

u  c 

c 

U  •-* 

W 

%J  u 

m 

44 

0  *J 

o 

0-0 

o 

*#4 

O' 

o 

c 

t A  C 

*. 

3 

o 

o 

0 

x: 

<* 

O' 

44 

-  x: 

in 

o 

« 

0 

O'  *J 

O'  V 

•*-0 

o> 

o 

a 

C  O' 

e  u 

\0 

•ft  1 

v  c 

•**  C 

0-0 

o 

© 

U  #J 

■X  0 

o 

% 

•ft 

• 

*4  -4 

O  3 

O' 

o 

© 

tt 

w 

«  Cf 

\ 

O' 

O' 

C 

rH 

0 

o 

* 

0 

U 

M»  U) 

V— 

o 

o 

fi 

0  0 

0— * 

■0 

04 

V  W 

o 

o 

O 

J2 

% 

0 

0 

U» 

o 

a 

a 

O' 

O' 

0 

•-4 

0 

o 

a 

0 

* 

• 

0 

X 

a 

o 

o 

4i 

C 

tel 

«—• 

o 

144 

\3 

M0 

0 

44 

0 

o  o 

© 

t*4 

0  V 

c 

XJ 

U  «-04 

ft** 

« 

fit.  -*  4J 

* 

1- 

0.  c 

© 

o 

o 

u 

0  0 

iQ 

r» 

0 

BO  U 

3  O  u 

w 

u 

#4  ft* 

o  5. 

0 

>  O'- 

> 

•>** 

« 

0 

«u 

0 

0  « 

0 

0 

<# 

a 

Q 

K 

•  ft** 

(Si 

N 

CM 

0 

©  •*• 

£  & 

« 

• 

•  i,  * 

* 

132 


\k 

'  % 

I  . 

|H 

]  k'  *’ 


W 

c  ui 

d,  «— 

E-  •-*> 
o 
c  ch 
*«•  \ 
m 


rH 

o 

VO 

Oi 

PM 

<N 

ro 

o 

o 

o 

o 

♦ 

• 

• 

• 

o 

o 

o 

o 

o 

eg 

o 

Ov 

eg 

m 

tr> 

VO 

rH 

in 

ro 

Lf> 

oooo  oooo 

oomji  o  in  r- 

lO  l£>  f"  f"*  rH  CTl  fM  •-< 

H  rt  H  M  n  N  «  « 


oooo  oooo 

tn  e  in  in  (OHO* 

r-  «*  in  p>-eii«#vO 

MriMft  in  m  m  m 


% 


a 


o 

C 

o 

o 

Q 

o 

o 

o 

• 

• 

» 

* 

CD 

CO 

03 

oeoo 

o  me  « 

a  c*  ©  a 

in  m  m  m 

in  «#  m  >n 

r-  r~  r*  !"* 

oooo 

OOOO* 

in 

in 

fM 

<M 

•4 

•4 

• 

« 

o 

o 

eov'iB 

QF>  O  f 

OOOO 

e  c?*  a  o 

o  o  o  o 

at  n  in 

•»«  —4  **4 

^ 

o 

o 

o 

«n 

rvi 

N 

• 

A 

o 

o 

134 


135 


w  U) 
c  •— 
a>  cm 
£-<  — ' 
o 
c  co 

+i 
\ 
r.  o 

4J  VO 
CP  +1 

c  \ 

4)  CM 


a 


x> 

« 

S- 


a 


4-1 

1 

o 

in 

in 

CM 

in 

a 

cn 

CO 

CM 

r- 

o 

CM 

> 

o 

cm 

m 

rH 

CM 

c 

0 

o 

r- 

CM 

CO 

<n 

00 

VO 

26 


o 

H 

O 

• 

o 


o 

o 


o 

H 

O 

« 

o 


VO 

c-t 

o 

4 

o 


in 

i—C 

o 

o 


•J  o 

CP 

OOOIA 

o  o  a  o 

o  o  o  o 

o  o  o  o 

o  o  a  o 

U 

vo  co  in  cm 

IflOOlD 

o  o  o  o 

Hinno 

o  m  in  o 

■H  M 

QU 

r-  <p  in  o 

CM  i-4  CO  CM 

in  co  aa  r» 

m  ro  h  o 

in  #-i  <o  cp 

4J 

H  rl  N  N 

CO  1*7  CO  CO 

CO  CO  CO  CO 

C-  CO  00  00 

vo  vo  n>  in 

w  e 
■-«  o 

u  •** 

4)  *J 

4J  «S 

U  44 
<0  C 

14  4) 

(0  >-< 

m 

in 

in 

in 

in 

XI  u 

CM 

CM 

CM1 

CM 

CM 

u  o 

X 

«-C 

*** 

1-4 

r4 

cH 

SN 

' 

• 

o 

• 

o 

4 

o 

• 

O 

• 

O 

o 

a 

o 

o 

O 

o 

o 

c> 

o 

o 

• 

4 

>* 

• 

4 

CO 

co 

co 

ao 

4»  iO£> 

COHOOO 

m  m  co  in 

4#!  CM  VP  © 

o  in  in  in 

©a  q  © 

co  co  «n  «# 

H  O  H  rt 

CP  SP  O'  O' 

p*  cn  vs>  vo 

in  in  m  in 

n>  p«  r»  **< 

o  o  ©  a 

O'  CMP  CP 

CP  O'  (P  CP 

ocdo 

o  o  o  o 

*•< 

OHHH 

^4 

in 

m 

c 

o 

© 

na 

CM 

tn 

in 

o 

H 

*M 

CM 

in 

• 

4- 

• 

4 

o 

O 

© 

O 

© 

Table  12  Characteristic  Length  in  Tension 
Ply  Orientation  [(0/90),) 
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APPENDIX  G 


Summary  of  Data  for  Loaded  Holes 

This  Appendix  contains  the  data  which  were  generated 
from  Fiberite  T30Q/1034-C  graphite  epoxy  laminates 
containing  loaded  holes.  The  Tables  in  this  Appendix  also 
contain  the  failure  strengths  and  failure  modes  calculated 
by  the  present  model  for  the  conditions  of  the  tests. 

Notations  used  in  Tables  15-29 


D  hole  diameter 

W  specimen  width 

L  specimen  length 

H  specimen  thickness 

E  edge  distance 

distance  between  two  parallel  holes 

Gy  distance  between  two  series  holes 

P  failure  load  under  tension 

?aVg  average  failure  load  under  tension 

Pc  calculated  failure  load 

H  experimental  failure  mode 

Mc  calculated  failure  mode 

T  tension  failure  mode 

B  bearing  failure  mode 

S  shearout  failure  mode 

T*  tearout  along  fiber  direction 

at  S45° 


(in) 

(in) 

( in) 
(in) 
(in) 

( in ) 
(in) 
(lbf) 
(lbf ) 
(lbf) 


Table  15  Data  and  Calculated  Values  for  Joints  Containing  a  single  Hole 

I ( 0/ ±45/90 ) _ ] 
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